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Abstract 

In this review, we give an overview of the experimental and theoretical advances 
in the physics of quantized vortices in dilute atomic-gas Bose-Einstein condensates 
in a trapping potential, especially focusing on experimental research activities and 
their theoretical interpretations. Making good use of the atom optical technique, 
the experiments have revealed many novel structural and dynamic properties of 
quantized vortices by directly visualizing vortex cores from an image of the density 
profiles. These results lead to a deep understanding of superfluid hydrodynamics of 
such systems. Typically, vortices are stabilized by a rotating potential created by a 
laser beam, magnetic field, and thermal gas. Finite size effects and inhomogeneity 
of the system, originating from the confinement by the trapping potential, yield 
unique vortex dynamics coupled with the collective excitations of the condensate. 
Measuring the frequencies of the collective modes is an accurate tool for clarifying 
the character of the vortex state. The topics included in this review are the mech- 
anism of vortex formation, equilibrium properties, and dynamics of a single vortex 
and those of a vortex lattice in a rapidly rotating condensate. 
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1. Introduction 



The achievement of Bose-Einstein condensation in trapped atomic gases at 
ultra-low temperatures has stimulated intense experimental and theoretical 
activity in modern physics, as seen by the award of the Nobel Prize in Physics 
in 2001 flCornelll 120021 iKetterle! l2002h . The Bose-Einstein condensate (BEC), 
a state of matter predicted by Einstein in 1925, is created by the condensation 
of a macroscopically large number of bosons into one of the eigenstates of the 
single-particle density matrix below the Bose-Einstein transition temperature. 
A remarkable consequence of the condensation is an extension of microscopic 
quantum phenomena into the macroscopic scale. This is an essential origin of 
superfluidity and superconductivity, in which macroscopically extended phase 
coherence allows a dissipationless current to flow. 

Superfluidity is closely related to the existence of quantized vortices. For 
weakly interacting BECs the superfluid velocity v is given by the gradient 
of the phase 9 of a "condensate wave function" v = (h/m)'V9 with the Planck 
constant H = h/2n and particle mass m. Since the wave function remains 
single-valued, the change in the phase around a closed contour must be an 
integer multiple of 2n. Thus, the circulation T around a closed contour is 
given by T — § v • dl — (h/m)q (q = 0, 1, 2, • • ■), which shows that circulation 
of a vortex is "quantized" in units of h/m. Realization of weakly-interacting 
atomic-gas BECs has provided an ideal testing ground to study the physics 
of quantized vortices; up to now, several experimental groups have reported 
many interesting results. This experimental work has been followed by con- 
siderable theoretical activity, leading to proposals and new problems to be 
tackled. (For a review of the early research stages of quantized vortices, see 
flFetterll2001a! u 

In this article, we review the physics of quantized vortices in atomic-gas BECs, 
especially focusing on the progress of the experimental research. We restrict 
ourselves to arguments on trapped condensates with inhomogeneous density 
profiles. The aim of this review is to stimulate further developments of this 
field. By reflecting on the history of the current research and on unresolved 
problems, we hope to encourage researchers in low temperature physics to 
investigate quantized vortices in this system. While quantiz ed vortices hay e 
been extensively studied in the field of superfluid helium ([Donnelly! If 9911 ). 
there has been a resurgence of interest in vortices in atomic BECs because 
of the following reasons. First, the diluteness of a gas yields a relatively large 
healing length that characterizes the vortex core size, thus enabling the visu- 
alization of vortex cores by imaging techniques characteristic of this system 
(see Sec. 2. 3). Because of this observational capability, the ability to manip- 
ulate a condensate wave function, and the tunability of the rotation over a 
wide range, these systems provide a unique approach to studying quantized 
vortices and their dynamics. Second, the finite size effect due to the trapping 
potential causes novel properties of vortices. Finally, multicomponent BECs 
provide new possibilities for studying unconventional vortex states that have 
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been studied in other fields of physics, such as superfluid 3 He, anisotr opic su- 
perco nductors, and theories in high-energy physics and cosmology (jVolovik 



2003|). 



Since this is the first time that the topic of atomic BECs has appeared in 
the Progress of Low Temperature Physics, we start with a basic introduction 
to ultra-cold atomic systems, including how a condensate is formed and how 
they are manipulated. Although it is desirable to refer to the experimental 
and theoretical studies of such systems in detail, we will only mention the 
basic ideas necessary to understand some important issues in experiment and 
theory because of space restrictions. These issues are described in Sec. 2. More 
detailed acc ounts can be found in the comprehensive t ext books by Pe thick 
and Smith (IPethickl 120021 ) . Pita evskii and S. Stringari (IPitaevskiil 120031 ) . and 
in the review paper by Leggett ( jLeggettll200ll ). In Sec. 3, we review the basic 
theory and experiments on quantized vortices in atomic BECs, addressing how 
vortices are created in this system and how they are detected. We also show 
the intrinsic mechanism of vortex nuclation and lattice formation in a trapped 
BEC. There are two interesting regimes classified by the rotation rate of the 
system: one at a slow rotation rate close to a critical rotation frequency where 
there is a single vortex, and another at high rotation rates for which a lattice 
of a large number of vortices is formed. The details of these two regimes are 
discussed in Sec. 4 and Sec. 5, respectively. Further interesting topics that 
cannot be explained sufficiently in this review and remaining future problems 
are presented in Sec. 6. We devote Sec. 7 to conclusions and outlook. 



2. Introduction to ultra-cold atomic-gas BECs 

While progress toward the achievement of Bose-Einstein condensation in a di- 
lute atomic gas had proceeded for the past few dec ades, researchers succeeded 



in creating a condensate in 1995 (j Anderson! 1 19951 ) . This realization brought 



great sensation in modern physics and opened a new research field combin- 
ing condensed matter physics and atomic, molecular, optical (AMO) physics. 
Here, we briefly summarize the basic introduction about a system of ultra- 
cold atomic-gas BECs, giving background information for understanding the 
experiments and theories of quantized vortices in the following sections. 

2.1. General information 

A typical system considered here is a collection of neutral atoms with particle 
number iV ~ 10 4 to 10 7 , trapped by a potential created by a magnetic field 
or an optical laser field. The density of the atomic gas is of the order of 
n ~ 10 14 cm -3 , which is lower than that of air on the earth (~ 10 19 cm -3 ). 
The transition temperature to Bose-Einstein condensation can be estimated 
from a dimensional analysis of the relevant physical quantities (m, n, H) as 
kT c ~ h 2 n 2 / 3 /m, which is in a range from 100 nK to a few /iK. At such low 
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temperatures, the gas phase cannot be a stable thermodynamic state and could 
in principle collapse to the solid phase. However, this relaxation is dominated 
by a three-body recombination process, which is a rare event for dilute and 
cold gases; the lifetime of the sample is thus long enough (of the order of a 
few seconds to a few minutes) to carry out experiments. 
In typical experiments, there are several steps toward the condensation of 
atoms. The first is laser cooling, achieved with three pairs of counter-propagating 
laser beams along three orthogonal axes. Subsequently, the precooled gas is 
confined in a trapping potential, described typically by a harmonic potential. 
In this stage, the temperature is of order 100 fiK, with 10 9 atoms. In the 
case of a trap created by a magnetic field, the atoms are trapped by the Zee- 
man interaction of the electron spin with an inhomogeneous magnetic field. 
Thus, atoms with electron spins parallel to the magnetic field are attracted to 
the minimum of the magnetic field (weak-field seeking state), while ones with 
electron spin antiparallel are repelled (strong-field seeking state). Laser cool- 
ing alone cannot produce sufficiently high densities and low temperatures for 
condensation. The second step, evaporative cooling (a process in some sense 
similar to blowing on coffee to cool it), allows the removal of more energetic 
atoms, thus further cooling the cloud. The evaporation is effected by apply- 
ing a radio-frequency magnetic field which flips the electron spin of the most 
energetic atoms. At the end of the process, the final temperature is about 100 
nK and about 10 4 -10 7 atoms remain. 

Experimentally, the atomic-gas systems are attractive, since they can be ma- 
nipulated by the use of lasers and magnetic fields. The cold gas is confined 
in a trap without microscopic roughness, as it is an extremely clean system. 
In addition, interactions between atoms may be affected either by using dif- 
ferent atomic species or by changing the strength of an applied magnetic or 
electric field for species that exhibit Feshbach resonance. A further advantage 
is that, because of low density, the microscopic length scales are so large that 
the structure of the condensate wave function may be investigated by optical 
methods. Finally, the mean collision time r co n ~ (anv ) _1 ~ 10~ 3 sec between 
atoms (a is the cross section and v the speed of atom) is comparable to the 
characteristic time of the collective mode, which prohibits a local equilibrium 
of the system. Thus, this system is ideal for studying nonequilibrium relaxation 
dynamics. 

2.2. Atomic Species 

The characteristics of BECs are mainly determined by atom-atom interac- 
tions, which depend crucially on the species of the condensed atoms. Most BEC 
experiments have been performed using alkali atoms because their ground 
state electronic structure is simple; all electrons except one occupy closed shells 
and the remaining electron is in an s orbital in a higher shell. This structure is 
well suited to laser-based manipulation because its optical transitions can be 
excited by available lasers and the internal energy-level structure is favorable 
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for cooling to very low temperatures. Since the first observation of BEC in 
atomic gases, BECs h ave been formed from nine differen t ele ments, includin g 
the alkali atoms 87 Rb flAndersonlll995h . 23 Na flDavislll995h. 7 Li flBradlevlll995h 



JljFriedlESh, 85 RblCornishl l200nh . 41 K f lModugnol boOlh . 133 Cs flWeber 



2001 



2003). and K (jRoatil 120071) . and the non-alka li atoms met astable He (Robert 



Dos SantoJT200lh . 174 Yb flTakasul l2003h . and 52 Cr flGriesmaierl l2005h 



87 Rb and 23 Na atoms are stable and have long lifetimes against inelastic col- 
lisional decay and are thus popular atomic species for BEC experiments. 



2.3. Detection 



Once a BEC has been created in a harmonic trap, it is probed for its prop- 
erties. This can be achieved either in situ, i.e., with the condensate inside 
the trap, or using a time-of-flight (TOF) technique. Al though in situ diag- 
nostics, such as nondestructive phase-contrast imaging (jAndrewd Il996l ). are 
valuable tools for some applications, the TOF technique is more often used in 
vortex experiments. The TOF technique involves switching off the trapping 
field (magnetic or optical) at time t — and taking an image of the BEC 
a few (typically 5 to 25) milliseconds later. Switching off the trap allows the 
sample to expand before applying the laser beam probe, because the probe 
is difficult to apply at high densities. Images of the sample are most often 
taken by absorption, i.e., shining a resonant laser beam into the atomic cloud 
and using a CCD camera to observe the shadow cast by the absorption of 
photons, from which can be determined the integrated atomic density. This 
method is inherently destructive since real absorption processes are involved 
by spontaneous radiation and the accompanying heating. 



2.4- Manipulation 



(i) Laser created potential 

Atoms in a laser field experience a force, due mainly to the interaction of the 
laser field with the electric dipole moment induced in the atoms. The force on 
atoms in a laser field is used in a variety of ways in BEC experiments. 
The character of the force is determined by the detuning given by A = 
^ias — (E e — E g ), where to>i as is the laser frequency and E g (E e ) the ground 
(excited) state energy of an atom. The force also depends on the laser-beam 
intensity 1$, given by Iq = eocT 2 /d 2 , where eo is the dielectric constant, c 
the speed of light, d an appropriately defined dipole matrix element for the 
transition in question, and T = h/r e with the lifetime of the excited state r e . 
In the limit r <C A, the change in energy of the atom in the laser field is 
Ai^ ascr (r) = (J(r)// )r 2 /A. A region of high laser intensity thus provides an 
attractive potential for A < ("red detuning") and a repulsive potential for 
A > ("blue detuning"). A red-detuned potential has been used as an opti- 
cal trap for atoms. A blue-detuned potential creates a potential barrier that 
separates a condensate and an impurity (obstacle) potential. The interference 
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pattern created by counter-propagating laser beams yields a periodic potential 
for atoms, called an optical lattice. 

(ii) Hyperfine state 

Atomic BECs can have internal degrees of freedom, attributed to the hyperfine 
spin of atoms. A hyperfine-Zeeman sublevel of an atom with total electronic 
angular momentum J and nuclear spin I may be labeled by the projection mp 
of total atomic spin F = I + J on the axis of the field B and by the value 
of total F, which can take a value from \I — J\ to 1 1 + J\. This is because 
the hyperfine coupling, which is proportional to I ■ J, is much larger than 
the typical temperature of an ultra-cold atomic system. The hyperfine state 
is denoted by \F, mp) with rap = —F, —F + 1, • • • , F — 1,F. The simulta- 
neous trapping of atoms with different hyperfine sublevels makes it possible 

ten c alled "spinor") BECs with internal degrees 
1998h. characterized by multiple order parame- 



to create m ulticomponent (o 
of fr eedom (Ho 19981. Ohmi 



ters flHalllll998l.lStengerl ll998 



Barrettll200ll . ISchmaliohannll2004 IChandl2004 



Kuwamotol 12004 ) . 

An external field can couple the internal sublevels of the atom and cause co- 
herent transition of the population. This coherent transition can be used to 
control the spatial variation of the condensate wave functions, resulting in an 
"imprinting" of a phase pattern onto the condensate. In most schemes, the 
spatial configuration of the field, the intensity and detuning of the laser fields, 
and the phase relationship between the different fields need to be carefully 
controlled to create the right phase pattern, that takes full advantage of the 
complex internal dynamics. 

(iii) Feshbach resonance 

A salient feature of cold atom systems is that fi eld-induced F eshbach resonance 
can tune the scattering length between atoms (Ilnouyeil200ll ) , which determines 
the atom-atom interaction. A Feshbach resonance occurs when a quasi-bound 
molecular state in a closed channel has energy equal to that of two collid- 
ing atoms in an open channel. Such resonances can greatly effect elastic and 
inelastic collisions such as dipolar relaxation and three-body recombination. 
Scattering near the resonance can be quantified by perturbation theory. To 
first order in the coupling between open and closed channels, the scattering 
is unaltered, because there are no continuum states in the closed channels. 
However, two particles in an open channel can scatter to an intermediate 
state in a closed channel, which subsequently decays to give two particles in 
an open channel. Considering such second-order processes, we can obtain the 
contribution to the scattering length as ~ (E op — i? cl ) _1 , where E op is the 
energy of the particles in the open channel and E c \ is the energy of a state in 
the closed channels. Consequently, there are large effects if the energy E op of 
the two particles in the entrance channel is close to the energy E c \ of a bound 
state in the closed channels. Therefore, coupling between the channels yields 
a repulsive interaction if the energy of the scattering particles is greater than 
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that of the bound state, and an attractive interaction if it is less. Since the 
energies of the states depend on external parameters such as the magnetic 
field, the resonances can be used to control the interaction between atoms. 

2. 5. Basic theory of trapped BECs 
2.5.1. The Gross-Pitaevskii equation 

From a theoretical point of view, a major advantage of weakly-interacting 
atomic-gas BECs is that almost all the atoms in the system occupy the same 
quantum state and the condensate may be described very well in terms of 
mean-field theory. This is in contrast to liquid 4 He, for which a mean-field 
approach is inapplicable due to the strong correlations induced by the inter- 
actions between the atoms. Bogoliubov's treatment of a uniform Bose gas at 
zero temperature provides a useful mean-field description of a condensate. Sub- 
sequently, Gross and Pitaevskii independently considered an inhomogeneous 
dilute Bose gas, generalizing Bogoliubov's approach to include nonuniform 
states, which includes quantized vortices. Such nonuniform states of a dilute 
Bose gas can be understood by considering the second-quantized many-body 
Hamiltonian 



H = j dr¥{r) + Kx(r) + \ J dv'¥ {v')V mt ^ ~ r'MO 



expressed in terms of Boson field operators ^(r) and ^(r) that obey Bose- 
Einstein commutation relations [^(r), &(r')] = 5(r-r'), [^(r),^(r')] = 
[yt(r), ^(r')] = 0. Here, V ex (r) is a trapping potential. The interparticle po- 
tential Vi n t is approximated by a short-range interaction V int ~ g5(r — r'), 
where g = 4iih 2 a/m is a coupling constant, characterized by the s-wave scat- 
tering length a, because only binary collisions at low energy are relevant in 
a dilute cold gas and these collisions are independent of the details of the 
two-body potential. 

In three dimentions, a remarkable feature of a dilute Bose gas at zero temper- 
ature is the existence of a macroscopic wave function \1/ (an "order parame- 
ter"). The macroscopic occupation of condensed particles makes it natural to 
write the field operator as a sum \P(r, t) = \l/(r, t) + 0(r, t) of a classical field 
\l/(r, t) that characterizes the condensate and a quantum field 0(r, t) repre- 
senting the remaining noncondensed particles. In order to derive the equation 
of motion for the order parameter, we write the time evolution of the operator 
^(r, t) = exp(iHt/H)ty(r) exp(— iHt/H) using the Heisenberg equation with 
the many-body Hamiltonian 



h 2 v 2 

+ V cx + g¥(r,t)V(r,t) 



2m 



*M). (2) 



To leading order, the Bogoliubov approximation neglects the noncondensed 
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contribution, giving the time-dependent Gross-Pitaevskii (GP) equation 



ih- 



Ot 



fi 2 X7 2 

— + Vi + S |*(r, ( )|* 



*(r,t) (3) 



for the condensate wave function \I/(r, t). The GP equation (3) can be used 
to explore the dynamic behavior of the system, characterized by variations 
of the order parameter over distances larger than the mean distance between 
atoms. This equation is valid when the s-wave scattering length is much smaller 
than the average distance between atoms, and the number of atoms in the 
condensate is much larger than unity. 

The ground state of a trapped BEC can be expressed within the formalism 
of the GP theory. We can write the condensate wave function as ^(r, t) = 
<l>(r)e -l/ ^ //? \ where 3>(r) obeys the time-independent GP equation 



h 2 V 2 

V ex + #|$(r) 



2m 



$(r)=/i$(r); (4) 



$ is normalized to the number of condensed particles J c?r|$(r)| 2 = iVo, which 
determines the chemical potenteial \i. Typically, studies of trapped atomic 
gases involve the dilute limit (the gas parameter n|a| 3 is typically less than 
10~ 3 , where n is the average density of the gas), so that depletion of the 



condensate is small with N' = N — N oc ^n\a\ 3 N <C N. Hence, most of the 
particles remain in the condensate such that iVo ~ N. The time-independent 
GP equation (4) is also derived by minimizing the GP energy functional 

E[V, **] = J dr** + V ex + ||^| 2 ) * = E kin + E tI + E- mt , (5) 

subject to the constraint of a fixed particle number N. This constraint is taken 
into account by the Lagrange multiplier method; we write the minimization 
procedure as 5{E — fiN)/S^* = 0, where the chemical potential \x is the 
Lagrange multiplier that ensures a fixed N. 

Equation (4) provides a starting point for studying the structure of a conden- 
sate in a harmonic confining potential V^ x = m(uj 2 x 2 + u 2 y 2 + u 2 z 2 )/2. This 

introduces the length scale Oho = Jh/mu with u = {uj x uj y uj z ) 1 ^ . Although the 



exact ground state can be obtained only by solving Eq. (4) numerically, an 
approximate analytic solu tion can be g ained when the interaction energy E- irsk 



is much larger than E kin ( ]Baymlll995l ). To see this argument, let us neglect 



the anisotropy of the harmonic potential and assume that the cloud occupies 
a region of radius ~ R, so that n ~ N/R 3 . Then, the scale of the harmonic 
oscillator energy per particle is ~ muj 2 R 2 /2 while each particle experiences an 
interaction with the other particles of energy ~ gN / R? . By comparing these 
energies, the radius is found to be R ~ ah (87rA^a/a ho ) 1//5 . The kinetic energy 
is of order h 2 /2mR 2 , so that the ratio of the kinetic to interaction (or trap) 
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energies is ~ (Na/ah )~^ 5 ■ In the limit Na/a ho ^> 1, which is relevant to 
current experiments on trapped BECs, the repulsive interactions significantly 
expand the condensate, so that the kinetic energy associated with the density 
variation becomes negligible compared to the trap and interaction energies. 
As a result, the kinetic-energy operator can be omitted in Eq. (4), giving the 
Thomas-Fermi (TF) parabolic profile for the ground-state density 

n(r) ~ |* TF (r)| 2 = - V ex (r)} , (6) 

where Q(x) denotes the step function. The resulting ellipsoidal density in 
three-dimensional (3D) space is characterized by two types of parameters: 
the central density n Q = fi/g and the three condensate radii i?| = l^jmuP- 
(j = x,y,z). The chemical potential /i is determined by the normalization 
/drra(r) = N as fi = (^/2)(15iVa/a ho ) 2/5 . 



2.5.2. The Bogoliubov-de Gennes equation 

The spectrum of elementary excitations of a condensate is an essential ingre- 
dient in calculations of the thermodynamic properties. To study the low-lying 
collective-excitation spectrum of trapped BECs, the Bogoliubov-de Gennes 
(BdG) equation coupled with the GP equation is a useful formalism. 
Let us consider the equation of motion for a small perturbation around the 
stationary state $, which is a solution of Eq. (4). The wave function takes the 
form \I/(r, t) = [<E>(r) + Uj(r)e~ tuJjt — v*(r)e lu,jt ]e~ llJ,t . By inserting this ansatz 
into Eq. (3) and retaining terms up to first order in u and v, we obtain the 
BdG equation: 



' £(r) -^(r) 2 











J 









(7) 



where C(r) = — ^ 2 V 2 /2m + V cx — fi + 2g\§(r)\ 2 , and ouj are the eigenfrequencies 
related to the quasiparticle normal-mode functions Uj(r) and Vj(r). The mode 
functions are subject to the orthogonality and symmetry relations / dr[uiU* — 
ViV*} = 5ij and / dv[uiV* — ViU*] = 0. 

Since the energy hujj of this quasiparticle is defined with respect to the con- 
densate energy, in Eq. (5) with the stationary solution $, the presence of 
quasiparticles with negative frequencies implies an energetic (thermodynamic) 
instability for the solution $. If there is energy dissipation in the system, the 
excitation of negative-energy modes lowers the total energy of the system and 
<3> relaxes to a more stable solution. We note that, since the matrix element of 
Eq. (7) is non-hermitian, the eigenfrequencies can be complex-valued. When 
there are complex- valued frequencies, small-amplitude fluctuations of the cor- 
responding eigenmodes grow exponentially during the energy-conserving time 
development. This is known as dynamical instability and is a main origin of the 
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generation of nonlinear excitations, such as vortices or solitons, in ultracold 
condensates. 



3. Vortex formation in atomic BECs 



In this section, we describe some basic properties of quantized vortices in 
trapped BECs and experimental procedures to create them. We also describe 
the nucleation mechanisms of quantized vortices in trapped BECs, including 
controlled schemes of phase engineering. Typically, experiments have used a 
smooth rotating potential, which is created by a laser or magnetic field, with 
a small transverse anisotropy to rotate the condensate. This potential induces 
a low-energy collective oscillation or shape deformation of the condensate. 
Such global motions of the condensate are responsible for the instability of 
the vortex nucleation, producing interesting nonequilibrium dynamics in the 
system. This is contrary to superfluid helium systems, where vortex nucleation 
occurs locally through roughness or impurities in the rotating container. 



3.1. Theoretical background 



As a simple example, let us consider the structure of a single vortex in a con- 
densate trapped by an axisymmetric harmonic potential V ex (r, z) = muj\{r 2 + 
\ 2 z 2 )/2 with aspect ratio A = oj z /uj^_. The condensate wave function with a 
quantized vortex line located along the z-axis takes the form $(r) = 0(r, z)e tq6 
with winding number q and cylindrical coordinate (r, 6, z). is a real function 
related to the condensate density as n(r,z) = <fi 2 . The velocity field around 
the vortex line is v s = (qh/mr)9. Equation (4) becomes 




d 2 | 1 d | d 2 \ | q 2 h 2 { . , 



= pl>. (8) 



Here, the centrifugal term q 2 h 2 /2mr 2 arises from the azimuthal motion of the 
condensate. Equation (8), solved numerically, gives the structure of the vortex. 
In the TF limit Na/a^o ^> 1, we can omit the terms involving derivatives with 
respect to r and z in Eq. (8), and the density can be obtained approximately 

as 



n(r, z) = n 1 - ~r - - q 2 ^ 6 1 - —r - — - q 2 — , (9) 




where n = \ij g is the density at the center of the vortex-free TF profile. We 
define the TF radius R\ = 2fj,/mw 2 _ and R 2 = 2fj,/mu 2 and the healing length 
£ = (h 2 /2mgno) l l 2 = (87ran ) _1 / 2 . Equation (9) shows that the condensate 
density vanishes at the center, out to a distance of order £, due to the cen- 
trifugal term £ 2 /r 2 (numerical solution shows that the density grows as r 2 
away from the center), whereas the density in the outer region has the form 
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of an upward-oriented parabola. Hence, the healing length £ characterizes the 
vortex core size; for typical BEC parameters, £ ~ 0.2 /im. In the TF limit, 
the core size is very small because £/R± = Hu±/2/j l = (IhNa/ 'ctho)~ 2//5 "C 1. 
Increasing the winding number q widens the core radius due to the centrifugal 
effects. 

The energy associated with a single vortex line is an important quantity to 
determine the stability of the vortex state. The dominant contribution to this 
energy is the kinetic energy of the superfluid flow by a vortex. The energy is 
estimated as 

Ei = / -mnv 2 dr ~ R z / vhnrdr = q 2 R z In -F , 10 

J 2 2j m \ s / 

where we assume the spatially uniform density n and the size along the z- 
axis R z . Since -Ei oc g 2 , vortices with q > 1 are energetically unfavorable. For 
example, the energy cost to create one q = 2 vortex is higher than that to 
create two q = 1 vortices. Therefore, a stable quantized vortex usually has 
q = 1, except for that in non-simp ly connected geometry, and we will mainly 
concentrate on the q = 1 vortex in the following discussions. In atomic BECs, 
however, such a multiply quantized vortex can be created experimentally by 



using topological phase imprinting (jShinll2004l ) and exhibits interesting disin- 



tegration dynamics, as discussed in Sec. 4.2.3. 

It is necessary to ensure the stability of a vortex in a trapped BECs against 
a non-vortex state to investigate its behavior. Imposing rotation on the sys- 
tem is a direct way to achieve stabilization. If the system is under rotation, 
it is convenient to consider the corresponding rotating frame; for a rotation 
frequency ft = flz, the integrand of the GP energy functional (5) acquires an 
additional term, 

e' = J dvm* {^~- + v cx + 9 -m 2 - vi/, (ii) 

where L z = —ih(xd y — yd x ). The corresponding GP equation becomes 



ih- 



dt 



h 2 X7 2 

-^— + V cx + g\^(r,t)\ 2 -nL 2 
2m 



*(r,t). (12) 



If there is a quantized vortex along the trap axis, (L z ) = Nh, so that the 
corresponding energy of the system in the rotating frame is E[ = E\ — NhQ. 
The difference between E[ and the vortex-free energy E' gives the favorable 
condition for a vortex to enter the condensate. Since E' is equal to the energy 
Eq in the laboratory frame, the difference is given by AE' — E[ — E' Q — 
Ei — Eq — NhQ. Thus, the critical rotation frequency Q c for the existence 
of an energetically stable vortex line is given by Q c = (E\ — Eq)/NH. Above 
the critical rotation frequency Q c the single vortex state is ensured to be 
thermodynamically stable. 
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To calculate E\ more quantitatively, it is neces sary to take into account the 
inhomogeneous effect of the condensate density ( Lundh 19971 ). In the TF limit, 
for a condensate in a cylindrical trap o> 2 = (an effective 2D condensate), the 
critical frequency is given by Q c = (2h/mR 2 [_) ln(0.888-R^/£). For an axisym- 
metric trap V ex (r, z), the critical frequency is fl c = (5h~/2mR 2 L ) ln(0.671-R^/£). 
fl c for a nonaxisymmetric trap is slightly modified by a sm all numerical fac- 
tor, and has been discussed analytically ( Svidzinsky 2000al ) and numerically 
flFederill999al ). 

When the rotation frequency is significantly higher than Q c , further vortices 
will appear in the form of a triangular lattice, in analogy to the Abrikosov 
lattice of magnetic fluxes in type-II superconductors. The nature of the equi- 
librium state changes, first to a state with two vortices rotating around each 
other, then to th ree vortices i n a triangle form, and subsequently to arrays 
of more vortices (iButtd Il999l . iFederl l2001bl ). The detail of this state will be 
described in Sec. 5. 



3.2. Vortex formation in a stirred condensate 



Rotation effects of atomic BECs were first studied based on the knowledge 
of deformed atomic-nuclei systems. A slight rotation of a deformed trap ex- 
cites so-called "scissors modes", which are closel y related to the irrotational- 
ity V x v s = of the superfluid hydrodynamics (lGuery-Odelmlll999l . iMarago 
20001 ). Although these experiments demonstrated the superfluidity of atomic 
BECs, more intuitive evidence can be gained by observing quantized vor- 
tices. However, since the first experimental realization of BECs, there have 
been technical hurdles to rotating a system. For atomic gases trapped by an 
impurity-free external potential, it was supposed that rotation of the potential 
could not transfer a sufficient angular momentum into the condensate. Thus, 
the first experimenta l detection of a vortex in an atomic BEC, by Matthews et 
al. ( jMatthewslll999al ) . was made by using a complicated ph ase imprinting tech- 
nique proposed by Williams and Holland (|Williamslll999l ). Subsequently, vor- 
tices have been created by stirring a condensate mechanically with an "optical 
spoon"; the firs t success of thi s "rotating bucket" experiment was reported by 



Madison et al. ( iMadisonl 120001 ). In the following, we detail the rotating bucket 
experiments. 



3.2.1. Experimental scheme to rotate condensates 

Madison et al. at Ecole Normale Superieure (ENS) succeede d in observing a 
vortex and a vortex lattice in an atomic BEC flMadisonll200oh . usin g a method 
simil ar to the observation of superfluid helium in a rotating bucket (jYarmchuk 
19791 ). A schematic illustration of their experimental setup is shown in Fig. 1. 
Since the magnetic trap is axially symmetric, its rotation cannot impart an 
angular momentum to the condensate. In this scheme, a laser beam is propa- 
gated along the z-axis of the condensate. This beam rapidly oscillates around 
the z-axis with a frequency much larger than the typical trapping frequency, 
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Fig. 1. Schematic illustration of the experimental setup used by the ENS group. 
A cigar-shaped condensate is rotated by an effective dipole potential made by the 
laser beam. The laser beam oscillates rapidly with a very large frequency around 
the z-axis, with an amplitude ~ 16 /im. The beam width is about 20 /mi. The dipole 
potential rotates with an angular frequency f2. 

which is effectively regarded to be as if two laser beams are located at an 
equilibrium position. Thus, the two laser beams break the axisymmetry of 
the condensate, allowing an angular-momentum transfer into the condensate. 
An optical spoon is then realized by rotating the two beams around the z- 
axis. Since the beam width is larger than the radial size of the condensate, 
the condensate is trapped in a trap combining an axisymmetric harmonic 
potential and a nonaxisymmetric harmonic potential created by a stirring 
laser beam. In a rotating frame, the combined potential can be written as 
(l/2)rW[[(l + e)X 2 + (1 - e)Y 2 } + (l/2)mu 2 z z 2 , where X and Y are the co- 



ordinates in the rotating frame, e 
parameter, and u± = 



X 



UJ 



Y 



)/(". 



\ + ujy) is the anisotropic 



uj 2 x + Uy)/2. By rotating this potential at a frequency 
Q, a vortex is formed above a certain critical value of Q after the equilibration. 
When Q is increased further, multiple vortices appear, forming a triangular 
lattice. The quantized vortices can be directly visualized as "dips" in the 
transverse density profile of the TOF image. 



Fig. 2. Typical density profiles of a rotating condensate taken by TOF measurement. 
The condensates contain approximately 16, 32, 80, and 130 vortices from left to 
right. From J.R. Abo-Shaeer et al., SCIENCE 292, 476 (2001). Reprinted with 
permission from A A AS. 

Following the experiments of the ENS group, other groups have also observed 
quantized vortices using slightly different methods under the concept of the 
rotating bucket. Abo-Shaeer et al. at Massachusetts Institute of Technology 
(MIT) observed a vortex lattice consisting of up to 100 vortices in a 23 Na con- 
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densate, as s hown in Fig . 2. 2 ^ Na condensates can be much larger than 87 Rb 
condensates (jAbo-Shaeerll200ll ) . Hodby et al. created a vortex lattice by rotat- 
ing the anisotropic magnetic trap directly without using a laser beam, which is 
similar to the rotating bucket experiment (lHodbyil2002h . This method has the 
advantage that a wider range of the anisotropic parameter e can be selected 
than for the optical spoon. Rotating an optical spoon made by multiple-spot 
l aser beams o r narrow focusing beams has also used for nucleating vortices 
flRamanllioOlh . 

In the above methods, the condensate was rotated by an external potential. 
Haljan et al. at Joint Institute for Laboratory Astrophysics (JILA), in con- 
trast, created a vortex stat e by cooling an initially rotating thermal gas in 
a static confining potential (IHaljanl l200ll ) . Thermal gas above the transition 
temperature was rotated by a slightly anisotropic trap. After recovering the 
anisotropy of the potential, the rotating thermal gas was evaporatively cooled 
until most of the atoms were condensed. Although the atom number decreased 
through the evaporative cooling, the condensate continued to rotate because 
the angular momentum per atom did not change and thus the vortex lattice 
was created. This method allows the investigation of the intrinsic mechanism 
of vortex nucleation, which is independent of the character of the stirring po- 
tential. In addition, since atoms can be selectively removed during the evapo- 
ration, spinning up of the condensate can be efficiently achieved by removing 
atoms extending in the axial direction, as opposed to the radial direction, and 
hence a BEC with a high rotation rate can be obtained. 



3.2.2. Theory of vortex nucleation and lattice formation 
(i) Surface mode instability 

The critical rotation frequency fl c indicates the energetic stability of the cen- 
tral vortex state and provides a lower bound for the critical frequency. Vortex 
nucleation of a non-rotating condensate occurs when the trap is rotated at 
a higher frequency than Q c , to overcome the energy b arrier that stop s the 
transition from the nonvortex state to the vortex state ( Isoshimalll999l ). The 
threshold of the rotation frequency for instability, leading to vortex nucle- 
ation, i s relat e d to the excitation of surface modes of the trapped conden sate 
flFederi Il999b1 . iDalfovol lioool . lAnriinl lioOll . IWilliamsl l200l ISm^uTal l2002ri ). It 
has been shown that, according to the Landau criterion for rotationally in- 
variant systems, the critical frequency is given by Q v = min(a^/£), where ue 
is the frequency of a surface mode with multipolarity I. Above Q v , negative- 
energy surface modes ap pear with high multipolarities in the spectrum of a 
non- rotating condensate (jlsoshimal Il999l . IDalfovol |2000| ). which may lead to 
vortex generation. The negative-energy mod es can grow on ly in the presence 
of dissipation, caused by e.g., thermal atoms (IWilliamsll2002l) . This mechanism 
occurred in the experiment of the JILA group (iHalj an! 120011 ). In this experi- 
ment, vortices were formed by cooling a rotating thermal cloud to below T c , 
where the surface modes were excited by the "wind" of the rotating thermal 
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cloud. 
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Fig. 3. Dependence of the vortex number (ang ular momentum p er atom) on the 
rotation frequency, (a) Result of the ENS group (|Chevvll2OO0l . l2OOlh . The parameter 
values are N = 2.5 x 10 5 , uj± = 2ir x 172 Hz, and an anisotropic parameter e (black 
dots: e = 0.01, white square: e = 0.02). The angular momentum was measured by 
the surface-wave spectroscopic technique described in Sec. 4.2.2. With increasing 
fi, the value of the angular momentum has a peak at O ~ uj±/\^2 = 122 Hz, 
a resonance frequency of the quadrupole mode. The width of this peak structure 
increas es as the trap ping potential becomes more anisotropic, (b) Result of the MIT 
group ( Ramanll200ll ). The parameter values are N ~ 10 and u>± = 2ir x 86 Hz. The 
arrows below the graph show the positions of the surface mode resonance u±/y/I. 
The i nset shows 2 -, 3-, and 4-point p otentials produced by a laser beam. (Taken 
from (|Chevvll200ih and (|E,amanll200ll ). Reprinted with permission from American 
Physical Society (APS).) 



The nucleation frequency fl v is insu fficient to expla i n the results of the groups 



using external stirring potentials (iMadisonl 120001 . lAbo-Shaeerl |2001| . iHodby 
20021 ). For example, in the case of the ENS group, the number of nucleated 
vortices has a peak near Q = 0.7u±, as shown in Fig. 3(a). These experiments 
confirm that instability occurs when a particular surface mode is resonantly 
excited by a deformed rotating potential. The optical spoon of the ENS group 
mainly excites the surface mode with I = 2 (quadrupole mode). In a rotating 
frame with frequency Q, the frequency of the surface mode is increased by — £Q. 
This resonance thereby occurs close to the rotation frequency O = ionjl. In the 
TF limit, the d ispersion relation for the surface mode is given by cue = \ftuj± 
( 1Stringarilll996l ). Hence, it is expected that the quadrupole mode with £ = 2 is 
resonantly excited at f2 = u±/^/2 ~ 0.707u;^. A theoretical study has revealed 
that, when the quadrupole mode is resonantly excited, an imaginary c ompo- 
nent appears in frequencies of fluctuations with high multipolarities (ISinha 
20011 ) . This indicates that dynamic instability can trigger vortex nucleation 
even at zero tempe rature. This sc heme is supported further by the experiment 
of the MIT group (jRamanll200ll ). where surface modes with higher multipo- 
larities (£ = 3, 4) were resonantly excited using multiple laser-beam spots, 
where the largest number of vortices were generated at frequencies close to 
the expected values £1 = u>±/VI, as seen in Fig. 3(b). 



3 



17 



Another mechanism to produce vortex nucleation is formulated by consider- 
ing elliptically-deformed stationary states of a BEC in a rotating elliptical 
trap. The stationary solutions can be obtained in t he TF limit by solving 
the superfluid hydrodynamic equations (jRecatil 120011 ). At low rotation rates, 
only one stationary solution is found ("normal branch"). However, at higher 
rotation rates, f2 > lu±/^ specifically, a bifurcation of the solutions occurs 
and up to three stationary solutions appear, with two solutions in an "over- 
critical branch" and one normal-branch solution. Above Q of the bifurcation 
point, one or more of the soluti ons become dynamically unstable for fluctua- 
tions with hig h multipolarities (ISinhall200ll ). which leads to vortex formation 
(jParkerl 120061 ) . Madison et al. followed these stationary states experimentally 
by adiabatically introducing trap ellipticity and rotation. T hey observed vo r- 
tex nucleation in the expected dynamically unstable region flMadisonll20"0ll ). 
An overview on the prob lem of vortex nucleation in a trapped BEC can be 
seen in the review paper (jGhoshl 12004 ). 



(ii) Formation dynamics of a vortex lattice 
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Fig. 4. (Left) Measurement of the time dependence of a (see text) when the stirring 
anisotropy is turned on rapidly from e = to e = 0.025 in 20 msec and Q = 0.7, 
and held constant for 300 ms. Five images taken at intervals of 150 ms show the 
transverse profile of the elliptic state and reve al the nucleation and ordering of 
the vortex lattice. (Taken from ( Madison 200ll ). Reprinted with permission from 
APS.) (Right) Comparing numerical simulation of the 3D GP equation with phe- 
nomenological dissipation (7 = 0.03) and the experimentally relevant parameters. 
The density profile in the bottom panels is integrated along the z-axis. 



Madison et al. directly observed nonlinear processes s uch as yortex nucle- 
ation and lattice formation in a rotating condensate (IMadisonl 1200 ll ). The 
left panel of Fig. 4 depicts the time development of the condensate elliptic- 
ity a = Q(R X — Ry)/(Rx + Ry), where Rx,y is the TF radius measured 
from the image. By suddenly turning on the rotation of the potential, the 
initially axisymmetric condensate undergoes a collective quadrupole oscilla- 
tion in which the condensate deforms elliptically. This oscillation continues 
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for a few hundred milliseconds with gradually decreasing amplitude. Then, 
the axial symmetry of the condensate suddenly recovers and concurrently the 
vortices enter the condensate from its surface, eventually settling into a lattice 
configuration. 

This ob servation has been well reproduce d by a sim ulation of the GP e quation 
in 2D (ITsubotal 12002k iKasamatsul 120031 ) and 3D (iKasamatsul l2005al ) space. 
The results shown in the right panel of Fig. 4 were obtained by a 3D sim- 
ulation with the parameters of the ENS group experiment. The simulation 
result shows that after a few hundred milliseconds, the boundary surface of 
the condensate becomes unstable and generates ripples that propagate along 
the surface, identified as invisible "ghost" vortices in the low-density surface 
region. The ripples develop into vortex cores, which enter the condensate. In 
these simulations, a dissipation term was introduced phenomenologically by 
rewriting the time derivative term of the GP equation (12) as (i — rfd/dt, 
which caused the lattice configuration to settle. Other works have simulated 
vortex lattice formation using dissipati on derived from the microscopic ap- 
proach, suc h as quantu m kinetic theory (IPenckwittl 120021 ) or the classical field 
formalism (lLobo 2004 ) . Long numerical propagation of the energy-conserving 
GP equatio n can cause c rystallization of a lattice through the vortex-phonon 
interaction (jParkerll2005l ) 



(iii) Vortex nucleation by a moving object 

Vortices can also be nucleated in BECs by a moving localized potential. Nu- 
merical simulations of the GP equation for a 2D uniform condensate flow 
around a circular hard-walled potential show that vo rtex-antivor tex pairs nu- 



cleate when the flow velocity exceed a critical value (lFrischlll992l ). In trapped 



BECs, a similar situation can be realized experime ntally using a narrow blue- 



detun ed laser potential, being studied theoretically (I Jacksonlll998l . ICrescimanno 



2000l ). In experiments by the MIT group, a repulsive laser beam was oscillated 



' Raman 


1999, 


Onofrio 


2000) 



the measurement of condensate heating a nd drag above a critical velocity was 
consistent with the nucleation of vortices (Jackson 2000l ). Focused laser beams 



moving in a circular p ath around the trap center can also stir t he condensate 



by nucleating vortices (ICaradoc-Davieslll999l. 12 000. Lu ndhll2003l ). This scheme 



was used in the experiment detailed in Ref. (jRamanll200ll ). where vortices were 
generated at lower stirring frequencies than the critical value given by surface 
mode instability. 



3.3. Phase engineering 



Atom optics techniques allow the controlled creation of vortices by imprinting 
a spatial phase pattern of the condensate wave function. Several ideas for the 
creation of vortices by this technique have been proposed theoretically, based 
on the coherent control of the time evolution of the wave function, instead of 



3 



19 



mech anical rotation f|Marzlinlll997l.ll998l.[Dun]ll998l. IWilliamdll999l iNakahara 



2000 



2003|) 



Andrelczykl l200ll . iDamski 



2001 



Nandi 



2004 iKapald l2005l . iMottonen 



3.3.1. Phase-imprinting method 

The first observation of a quantized vortex in an atomic-gas BEC was achieved 
in a two-component BEC consisting of 87 Rb atoms with hyperfine spi n states 
\F = l,m F = -1) = |1) and \F = 2,m F = 1) = |2) flMatthewslll999ah . which 
were confined simultaneously in almost identical magnetic potentials. Since 
the scattering lengths between atoms of |1) and |1), |2) and |2), and |1) and 
1 2) are all different, the two states are not equivalent, and the two-component 
condensate is characterized by two-component order parameters. 
In this experiment, condensed atoms are initially trapped in one state, say, 
the |1) state. Then, a two-photon microwave field is applied, inducing coherent 
Rabi transitions of atomic populations between the |1) state and the |2) state. 
For a homogenous system in which both components have uniform phases, 
interconversion takes place at the same rate everywhere. However, the time 
variation of the spatially inhomogeneous potential changes the nature of the 
interconversion. This is a key point of the phase- imprinting method for vortex 
creation. 

The underlying physics can be understood by considering a co-rotating frame 
with an off-centered perturbation potential at the rotation frequency Q'. In 
this frame, the energy of a vortex with one unit of angular momentum is 
shifted by hVL' relative to its value in the laboratory frame. When this energy 
shift is compensated for by the sum of the detuning energy of an applied mi- 
crowave field and the small chemical potential difference between the vortex 
and non-vortex states, a resonant transfer of population can occur. Experi- 
mentally, the rotating perturbation is created by a laser beam with a spatially 
inhomogeneous profile, rotating around the initial nonrotating component, say 
|1). By adjusting the detuning and Q', the |2) component is resonantly trans- 
ferred to a state with unit angular momen tum by precisel y controlling the 
time when the coupling drive is turned off ([Williams! Il999l ) . This procedure 
results in a "composite" vortex, where the |2) component has a vortex at the 
center, whereas the nonrotating |1) component occupi es the center and w orks 
as a pinning potential that stabilizes the vortex core (IKasamatsul l2005bl ) . 
As shown in Sec. 2.4, a far-off-resonant laser beam can create an external 
potential V[ as (r) in the condensate. By directly applying a laser pulse with 
an inhomogeneous intensity to the condensate, the condensate phase can be 
modulated. This can be easily understood by observing the evolution of the 
phase by inserting the form \I/(r, t) = |\I/(r, t)\ e^ r '*) into Eq. (3). When the 
laser intensity is much stronger than the other terms and the duration of the 
pulse r is sufficiently short, the evolution of the phase is governed by 8(r, t) = 
— ^ 1 Jo dtV\as(T, t). Since the potential amplitude of V[ as is proportional to 
the laser intensity, a suitable spatial variation of the intensity can imprint the 
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phase into the condensate (lAndrelczvkll2001l). This phase- imprinting te chnique 
has been used to create a dark soliton ( Burger 1999 . Denschlag 2000l ). which 
is a topological excitation with a complete density dip across which the phase 
changes by tt. It is known that a dark soliton in a dimensional space larger 
than 2D experiences dynamical instability, called "snake instability". This 
i nstability causes the decay of the dark solitons into a form of a vortex ring 
riAndersorJl200ll . button! l200lh . 



3.3.2. Topologic al vortex formatio n 



Leanhardt et al. (ILeanhardt 
phase imprinting" (INakahara 



2002, 12003T) used a method called "topological 



20001 . Ilsoshimall2000l . lQgawall2002h to create a 



vortex in a trapped BEC. In this experiment, Na condensates were prepared 
in either a lower, \F, thf) = |1, — 1), or upper, |2, +2), hyperfme state and 
confined in a Ioffe-Pritchard magnetic trap, described by B = B'(xx — yy) + 
B z z. A vortex was created by adiabatically inverting the axial bias field B z 
along the trap axis. 

To interpret the mechanism, consider an alkali atom with a hyperfine-spin 
\F\ = 1. The order parameter has three components ty±i and \l/o correspond- 
ing to F z = ±1, 0, respectively. The basis vectors in this representation are 
{|±), |0)}. We introduce another set of basis vectors \x), \y) and \z), which 



are defined by F x \x) 
previous vectors as 



F y \y) = F z \z) = 0. These vectors are related to the 



\z). When the 



±1) = T(l/V2)(\x)±t\y)) and |0> 
z-axis is taken parallel to the uniform magnetic field, the order parameter of 
the weak field seeking state takes the form = if) and — ^i — 0, or 
in vectorial form as SI/ = {ip/y/2) (x — iy). When the magnetic field points in 
the direction B = (sin (3 cos a, sin (3 sin a, cos j3), a rotational transformation 
with respect to the Euler angle (a, [3, 7) gives ^ = (■?/>/ v^2)e n (m + in), where 
m = (cos f3 cos a, cos (3 sin a, — sin (3) and n = (sin a, — cos a, 0). The unit vec- 
tor 1 = rh x n = (cos a sin (3, — sin a sin f3, — cos (3) is the direction of the spin 
polarization. The three real vectors {1, m, n} form a triad, analogous to the 
order parameter of the orbital part of superfluid 3 He. The same amplitudes in 
the basis {|0), |±)} are ^ = (^/2)(1 - cos f3)e~ ia+ ^ , # = -(V»/V% sin/?e i7 , 
and = (VV2)(1 + cos(3)e ta+t ^. 

When the field B z is strong compared to the quadrupole field, the trapped 
condensate has an order parameter \I> = (■?/;/v / 2)(x — iy) without vorticity. 
This configuration corresponds to /3 — and 7 = —a = <p, where <fi is the 
azimuthal angle. Then, B z is adiabatically decreased, where the adiabatic 
condition is required for atoms to remain in the weak field seeking state so that 
1 is always antiparallel to B. In the final step, the external field B z is gradually 
increased in the opposite (—z) direction. Then, the 1-vector points so that (3 = 
ii. Substituting these angles into \l/±i and we obtain = = and 
^1 = if)e 2l< ^, which corresponds to a vortex with win ding number q = 2. This 
resul t can be reinterpreted in terms of Berry's phase (jOgawall2002l . ILeanhardt 



20021 ) (which is why it is called "topological phase imprinting"). When the 
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hyperfine spin is F in general, we obtain a vortex with a winding numb er 2F 
since and ^ f hav e phases F(a + 7) and F(— a + 7), respectively (IShin 
200ilKumakurall2006h . 

When the bias field vanishes during the inversion (B z = 0), a spin texture 
known as cross disgyration appears in the Ioffe-Pritchard trap. Here, the angle 
(5 increases from to 7r/2 and 7 and —a are identified with 0, where the 
1-vector aligns with a hyperbolic distribution around the singularity at the 
center. This texture has a nonvanishing vorticity n when 7 = ncj). This spin 
texture has been observed as a coreless v ortex composed of three-component 
order parameters ^±1,0 of a spinor BEC (ILeanhardtl 120031 ) . 



3.3.3. Stimulated Raman process 

Some papers proposed generating vortices in a BEC using stimulated Raman 
processes with c onfigurations of opt i cal fields that have orbital angular m o- 
mentum (OAM) flMarzlinll 19971 . Il998l . |Pu"mlll998l . lNandll2004j . [Kapal3l2005h . A 
light beam with a phase singularity, such as a Laguerre- Gaussian (LG) beam, 
has a well-defined OAM along its propagation axis. The set of LG modes 



LGj(r, 



1 fy/2r 
ir(\l\ +p)w \ w 



2p 



2r2 



-r 2 1 'w 2 +U4> 



(13) 



defines a possible basis set to describe paraxial laser beams, where wq is the 
beam width, / the winding number, and p the number of radial nodes for radius 
r > 0. Each photon in the LGj, mode carries OAM lh along its direction of 
propagation. 



A group at NIST (1 Anderson! 12006k lRyull2007l ) used a 2-photon stimulated Ra- 
man process with a Gaussian laser beam propagating along +x and a LGq 
beam, carrying h of OAM, propagating along — x. Interference of counter- 
propagating Gaussian beams generates a moving sinusoidal optical dipole po- 
tential, which can give a directed linear momentum (LM) to Bose-condensed 
atoms via Bragg diffraction. The potential generated by interference of the 
counter propagating LGq and Gaussian beams is not sinusoidal but corkscrew- 
like, due to the radial intensity profile and the helical phase of the LGq beam. 
Diffraction off this corkscrew potential produces a vortex state with a center- 
of-mass motion, where atoms that absorb a photon from one beam and simu- 
latedly emit a photon into the other beam acquire both LM and OAM differ- 
ence of the beams, which in this case was 2hk (k the photon wavevector) and 
h, respectively. 

They generated vortices of higher charge by transferring to each atom t he an- 
gular momentum from several LGq photons ( Anderson! 120061 . lRyull2007l ). This 
experiment directly demonstrated that the OAM of a photon is transferred 
coherently to an atom in quantized units of h. In some situations it might be 
desirable to generate rotational states with no net LM. This could be accom- 
plished by using an initial Bragg diffraction pulse to put atoms in a non-zero 
LM state from which they could subsequently be transferred to a rotational 
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state with zero LM. 



4. A single vortex in an atomic BEC 



In this section, we concentrate on the problem of a single vortex state in a 
trapped BEC. As described in Sec. 3.1, vortex stability is ensured by the ro- 
tation of the system. Studying the motion of a vortex line is the first step 
towards understanding superfluid hydrodynamics in such a system. Trapped 
BECs are mesoscopic systems in the sense that the healing length is not signif- 
icantly smaller than the sample size. Thus, vortex dynamics have noticeable 
effects on the collective excitation of the condensate, in contrast to the case 
for traditional superfluid helium systems. 



4-1. Equilibrium properties 



The solution of the GP equation shows that a vortex has a core with a size 
of the order of the healing length £, in which the condensate density is zero. 
When Q = Qz, vortices are identified by a density dip in the transverse 
density distribu tion (in the xy plane). TOF observations by the ENS group 
(lMadisonll2000l ). however, show that the density is not completely zero in the 
dip. This result implies that the vortex line is not necessarily straight, because 
the condensate density along the z-axis (rotation axis) is integrated for the 
transverse image. Surprisingly, such vortex bending remains stationary in the 



ground state of a cigar-shaped condensate (lGarcia-Ripollll200ll . lAftalionll2002 



Modueno 


2003. 


Aftalion 


2003) 



Eviden ce of vortex bend ing in the ground state was observed by the ENS 
group (IRosenbuschl 120021 ) . They prepared a single vortex state slightly above 
Q c and equilibrated it for a sufficient long time. In the TOF measurements, 
two imaging beams were aligned along the y and z directions and probed the 
atom distribution simultaneously [Fig. 5 (a)]. The transverse image in Fig. 5 
(b) shows the vortex line, corresponding to the lower atom density; it is not 
straight and has the shape of a wide "U" . Figure 5 (c) shows the decay of a U 
vortex for which the angular momentum has decreased significantly compared 
to Fig. 5 (b). In the longitudinal view, we can see a vortex off-center. In the 
transverse view, a narrow U can be seen. 

This result is supported by theoretical analysis based on the 3D GP equa- 
tion with experimentally appropriate parameters (j Garcia- RipolfeOQll . lAftalionl 



20021 . lModugnoll2003l . lAftalionl 120031 ) . where the ground state with a rotation 
f2 = Qz was calculated by minimizing the energy functional in Eq. (11). 
The central vortex is generally bent if the trap aspect ratio A = uj z /uj± is 
much less than unity. A simple physical picture of the bending can be gained 
by viewing a ci gar-shaped condensate as a series of 2D sheets at various z 
(jModugnd 120031 ) . For each sheet, there is a corresponding 2D vortex stability 



problem with a critical frequency f2 



2D I 



(see Sec. 3.1) above which a cen- 
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Fig. 5. (a) Schematic of the imaging method of a vortex line. Two beams image the 
atom cloud simultaneously along the longitudinal (z) and transverse (y) directions 
of the initial cigar. In (b)-(d), the left column shows the "longitudinal" view along 
z, representing the atom distribution in the xy plane. The right column depicts the 
"transverse" view along the y direction, representing the atom distribution in the 
xz plane. The images were taken after equilibration times of (b) 4 s, (c) 7.5 s, (d) 
and 5 s. In the atom distribution of the transverse image, the ellipticity is inverted 
with respect to the initial cigar form , caused by transverse expansion during the 
TOF. (Taken from (|Rosenbuschll2002h . Reprinted with permission from APS.) 



tered vortex is the stable solution. Since the effective 2D chemical potential 
is H2d{z) = H — rnuj 2 z 2 /2, the radius of the 2D condensate at z becomes 
R±(z) = 2fi 2 D(z)/muj_ = i?j_(0) - \ 2 z 2 . Thus, tfPiz) oc R± is a decreasing 



function from z = to 



R z . For a given rotation frequency f2, the vortex 



line minimizing the total energy is well centered for \z\ < z c and is strongly 
bent for \z\ > z c , where Q 2D (z c ) = Q. This bending is a symmetry-breaking 
effect, which does not depend on the presence of rot ating anisotropy an d 
which occurs even in a completely axisymmetric system ( Garcia- Ripollll200ll ) . 
A precursor of this bending effect can be fou nd in the excitation spectrum o f 
a condensate with a centered straight vortex (jSvidzinskyll2000bl . lFederll2001al ). 
in which negative-energy modes localized at the core (so-called "anomalous 
modes") appear with increasing A. As these modes grow, the central vortex is 
pushed outward. This indicates that the bending instability needs a dissipation 
mechanism and that it takes a long time at low temperatures. 
Figure 5 (d) shows a vortex line in the shape of an "S" , which can be regarded 
as a U vortex with a half part rotated by 180°. An S single vortex can also be 
found by numerical simulation as the stationary sta te of an elongated conden- 
sate for a given rotation fr equency (lAftalioru 120031 ) . having an energy higher 
than that of the U vortex (jKominead 120051 ) . 
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4-2. Dynamical properties 



4-2.1. Precession and decay of an off-centered vortex 

Precession of a vortex core off-center in a condensate is a simple example of 
vortex motion. Core precession can be described in terms of a Magnus force 
effect. A net force on a quantized vortex core creates a pressure imbalance, 
resulting in core motion perpendicular to both the force and the vortex quan- 
tization axis. In the case of trapped BE Cs, these net forc es can be caused 
by either co ndensate density gradients (ISvidzinskvl l2000al Jbl. IJacksonl Il999 



McGedl200ll ) or drag due to thermal atoms (lFedichevlll999l ). The former may 



be thought of as a sort of effective buoyancy. Typically, the total buoyancy 
force is towards the condensate surface and the net effect is a precession of 
the core around the condensate axis via the Magnus effect. The latter causes 
radial drag and spiraling of the core towards the condensate surface due to 
energy dissipation and damping processes. 



Core precession has been investigated in detail by the JILA group ([Anderson 



2000h . Starting with a composite vortex created by the phase-engineering 
method of Sec. 3.3, they selectively removed components filling the vortex 
core with resonant light pressure. In the limit of complete removal, a single- 
component vortex state with a bare core can be obtained. The vortex was 
off-center because of the instability of the formation process. The precession 
frequency was determined from the vortex position taken directly from the 
density profile. The vortex core precessed in the same direction as the vortex 
fluid flow around the core. The obtained result of 1. 8 Hz agrees well w ith an- 



alytical results based on the Magn us force pictur e ( Svidzinsky 2000a bl) an d 
more precise numerical simulations (lJacksonlll999l . lMcGeell200ll . lFederll2001ah . 



In some results, the vortex core disappeared from the observed images dur- 
ing the time evolution. However, these were not associated with the decay 
of vortices because there was no evidence of radial spiral mo tion due to en- 
ergy dissipation, which could be caused by the t hermal drag (IFedichey I1999T) 
or the sound radiation from a moving vortex flLundhl |2000L IParkeri l2o"o4T ). 
Subsequently, an experiment using the surface-wave spectroscopic technique 
revealed that vortices were actually present for a long time in the condensate 
( lHaljanll200ll ). The main cause of the disappearance was the tilting motion of 
a vortex, the lowest odd-order nor mal mode of a sing le-vortex state, which is 
sensitive to small trap anisotropy (jSvidzinskyl l2000bl ) . 



4.2.2. Vortex dynamics coupled with collective modes 

Vortex dynamics are greatly affected by the overall collective motion of the 
condensate because of the mesoscopic nature of the system. Here, we detail 
several interesting results of coupled dynamics. The determination of the fre- 
quency of the collective modes allows precise measurement of the angular 
momentum carried by the quantized vortices. 
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(i) Transverse quadrupole mode 

The collective modes of a trapped condensate can be classified by expressing 
the density fluctuations in terms of polynomials of degree in the Cartesian co- 
ordinates (sci, £2,^3) = (x,y,z). The quadrupole modes are characterized by 
a density fluctuation with a polynomial of second order, e.g., 5n = J2Pij x i x j, 
which gives six normal modes. In an axisymmetric harmonic potential, linear 
combinations of the diagonal components p xx , p yy , and p zz describe three nor- 
mal modes: one transverse mode with m z = 2 and two radial-breathing modes 
with m z = 0, where m z is the projected angular momentum on the symmetry 
axis. The remaining three normal modes are associated with the off-diagonal 
components p xy , p yz , and p zx , which are scissors modes (IGuery-Odelinl Il999l . 
MaragdlioOOf ). 

The first study w as done for th e excitation of two transverse quadrupole modes 
with m z = ±2 (jChevyl 12000 ) . When a vortex is present in a condensate, 
the frequencies of the m z = ±2 quadrupole modes increase by ui+ — o>_ = 
2(L Z ) /mN(x 2 + y 2 ) because of broken rotational symmetry (jZambellil Il998h . 
where ( ) stands for the average within the condensate. The increase causes 
precession of the eigenaxes of the quadrupole mode at an angular frequency 
6 = (uj + — ujj)/2\m z \. By measuring the angular velocity of this precession, 
we can determine the mean angular momentum (L z ) of the condensate. This 
spectrosc opic method has also been used to characterize the tiltin g motion of 



a vor tex (jHaljan!l200ll ) and the winding number of a single vortex (ILeanhardt 
2002( 1 . 



Excitation of the transverse quadrupole mode yields further interesting vor- 
tex dynamics. The ENS group observed that when the superposition of the 
m z = ±2 quadrupole modes are excited with equal amplitudes, the osc illation 
of th e m z = —2 mode decays faster than that of the m z = +2 mode (IBretin 
20031 ) . A possible physical origin of this phenomenon is that the m z = —2 mode 
decays to Kelvin modes through a non-linear Beliaev pr ocess. This is sup - 
ported by theoretical analysis based on the BdG equation (jMizushimal 120031 1 . 
Kelvin modes correspond to long-wavelength helical traveling waves along a 
vortex line with a dispersion relation c^k — (frk 2 /2m) ln(l//c£) <C 1). Ac- 
cording to the Kelvin-Helmoltz theorem, the Kelvin modes rotate always in 
the sense opposite to the vortex velocity field. Consequently, the angular mo- 
mentum of a quantum of a Kelvin mode associated with a singly-quantized 
vortex is —h. Because of the negative angular momentum with respect to the 
vortex winding number, this mechanism is effective only for the m z = —2 
mode. From the energy uj_ 2 = 2co>k and angular momentum conservation, an 
excitation of the quadrupole mode m z = —2 can decay to a pair of Kelvin 
waves with wave vectors k and —k, while angular momentum conservation 
forbids the decay of the m z = +2 mode. 



(ii) Gyroscope motion 

What happens when the other quadrupole modes are excited in a condensate 
with a vortex? The Oxford group studied the response of a condensate with 
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a vortex when the xz or yz scissors modes are excited (IHodbyi 120031 ) . Similar 
to the case of transverse quadruople modes, in the presence of the vortex, the 
plane of oscillation of a scissors mode precesses slowly around the z axis. In 
polar coordinates, the scissors oscillation is in the 9 direction and the pre- 
cession is in the <fi direction, as shown schematically in F ig. 6. This can be 
regarded as a kind of gyroscope motion of the vortex line (jStringaril 120011 ). 




Fig. 6. Left: schematic picture of gyroscope motion of the experiment of Ref. (jHodbv 
2003). The scissors mode involves a fast oscillation of the small angle 9 between the 
condensate normal axis and the z axis. When a vortex is present, the plane of this 
oscillation (initially the xz plane with <f> = 0) slowly precesses through angle <j> about 
the z axis. Right: data of the evolution of the tilt angle projected onto the xz 
plane, when the scissor s mode is init ially excited (a) in the xz plane and (b) in the 
yz plane. (Taken from ( Hodbv 2003). Reprinted with permission from APS.) 



The relationship between the precession rate and (L z ) can be derived by con- 
sidering the scissors mode as an equal superposition of two counter-rotating 
m z = ±1 modes. These modes represent a condensate tilted by a small angle 
from the horizontal plane rotating around the z axis at the frequency of the 
scissors oscillation, uj± = u sc . The symmetry and degeneracy of these modes 
are also broken by the axial angular momentum (L z ). By applying a similar 
argument as that for transverse quadrupole modes f Stringaril 2001 ). the pre- 
cession rate is related to the frequency splitting, u + — U- = (L z )/mN(x 2 + z 2 ), 
allowing the angular momentum (L z ) to be determined. 
The precession associated with gyroscope motion can be seen in the results 
of Fig. 6(a) and (b), corresponding to a slowly varying oscillation component. 
The pattern of increase and decrease of the amplitude is exchanged between (a) 
and (b), with different directions of the excitation. This is clear evidence of a 
slow precession along the ^-direction. The results also show that the motion of 
the vortex core exactly follows the axis of the condensate. These observations 
can be well r eproduced by direct numerical simulations of the 3D GP equation 
(INilsenl 120031 ). From the precession rate, the measured angular momentum per 
particle associated with a vortex line was found to be 1.07 h ± 0.18/1. 



4-2.3. Splitting of a multiply quantized vortex 

The energy cost to create a q > 1 vortex is less favorable than q single- 
quantized vortices, as seen in Eq. (10). This raises an interesting question 
as to what happens when such an unstable vortex is created. As it happens, 
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topological phase imprinting, shown in Sec. 3.3.2, can be used to create an 
unstable vortex. 

The stability characteristics of a multiply quantized vortex in a trapped BEC 
exhibit an interesting interaction dependence because of the finite size effect 
( lPulll999l ). BdG analysis for a cylindrical system reveals complex eigenvalue 
modes, which implies that a multiply quantized vortex is dynamically unsta- 
ble. For a vortex with a winding number q, angular momentum conservation 
leads to constraints on the normal mode functions Uj(r) = Uj(r)e l( - K:i+q ^ e and 
Vj(r) = Vj(r)e t<yKj ~ qS>e , where Kj denotes the angular momentum quantum num- 
ber of the mode. For q = 2 there are alternating stable and unstable regions 
with respect to the interaction parameter an z = a J | 2 dxdy; the first and 
second regions appears for < an z < 3 and 11.4 < an z < 16. Numerical sim- 
ulations demonstrate that when a system is in an unstable region, a doub ly 
quantized vortex decays into two singly quantized vortices (iMottonenl 120031 ) . 
An experiment by the MIT group studied the splitting process of a d oubly 



quan tized vortex and its characteristic time scale as a function of an z=0 (I Shin 



20041 ). The results show that a doubly quantized vortex decays, but that the 



lifetime increases monotonically with an z=Q , showing no periodic behavior. 
This contrary to the above theoretical prediction. Recent numerical studies 
of 3D GP equations reveal this mysterious observation, emphasizing that the 
detailed dynamical behavior of a vortex along the entire z-axis is relevant for 
characterizing the splitting process in an elongated condensate (IHuhtamaki 



2006al . iMunoz Mated 120061 ) . 



The trigger for splitting instability is likely to be gravi tational sag during the 
formation process with a reversing axial bias field B z (IHuhtamakil l2006a[ ) . In 
experiments, absorption images were restricted to a 30-/im thick central slice of 
the condensate to increase the visibility of the vortex cores. The experimental 
results in Fig. 7(a) shows that the fastest decay occurs at an z ~ 1.5, consistent 
with the BdG and numerical analysis. As the particle number increases, the 
first instability region (0 < an z < 3) moves progressively away from the central 
slice toward the edges of the condensate because of the trapping potential. As 
a consequence, the splitting instability of the vortex core has to propagate 
from those regions to the center. This process is responsible for the monotonic 
increase in the lifetime for an z > 3. According to the theory, a second minimum 
is expected about an z= o ~ 13.75. Even though no such minimum occurs, 
a change in the slope of the predicted cu rve at an z =n — 13-75 is seen as 
an z= ft enters the second instability region (IHuhtamakil l2006al . IMunoz Mateo 



20061 ). Figures 7 (b)-(e) show the time evolution of the splitting process for 



an z= Q = 13.75. The first and second instability regions correspond to the 
shaded zones in Fig. 7 (b). This clearly shows that at t — 25 ms, the splitting 
process has already begun in both the edges and the center of the condensate. 
The different precession frequency along different z slices causes inter- winding 
of two single-quantized vortices. However, the two vortex cores near the central 
slice still overlap (Fig. 7 (c)) and thus cannot be experimentally resolved until 
much longer times. At t ~ 70 ms (Fig. 7 (d)), the vortex cores begin to 
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Fig. 7. Splitting process of a doubly quantized vortex, (a) Experimental data com- 
bined with theoretically predicted splitting times as a function of an z= o for a 4% 
quadrupolar perturbation acting during 0.3 ms. The splitting is identified by the 
number of visible vortex cores from the density profile of the axial absorption images 
taken from a 30-//m thick central slice of the condensate, (b)-(e) Time evolution 
of the splitting process for an z= o = 13.75, obtained by a 3D simulation of the GP 
equation. The shaded zones in (b) indicate the instability regions. The correspond- 
ing axial absorption images of the central slice are also sho wn at the bottom o f 
the figure. The lengths are in units of 6.05 fim. (Taken from ( Muhoz Mateo 20061 ). 
Reprinted with permission from APS.) 



disentangle so that they can be unambiguously resolved at t — 75 ms (Fig. 
7(e)). Thus, there is no contradiction with the theoretical prediction. 
The physical origin of the periodic appearance of the unstable region is anoma- 
lous modes with negative eigenvalues. When the eigenvalue of positive-energy 
modes coincides with the absolute value of the eigenvalue of negative-energy 
modes, a complex eigen value mode are produced through mutual annihilation 
of the se two excitations (jSkryabinll2000l . lKawaguchill2004l . I Jacksonll2005l . iLundh 
20061 ); the total angular momentum of these excitations is also vanished. Since 
the an z dependence of the negative-energy eigenvalues is very different from 
that of positive-energy ones, the above matching condition can be satisfied in 
sequence with increasing an z , which results in the periodic appearance of the 
complex-eigenvalue modes. Thus, splitting instability of a multiply quantized 
vortex can be suppressed for a particular trap asymmetry and interaction 
strength because the coll ective excitations d epend strongly on the character 
of the trapping potential (IHuhtamakil l2006bl ) . It has also been predicted that 
multiply quantized vor tices can be st abilized by introducing a suitable local- 
ized pinning potential flSimulall20 02a) or non-simply connected geometry such 
as quartic confinement (ILundhl 120021 ) . Very recent ly, splitting dy namics of a 
quadruply quantized (q = 4) vortex was reported (llsoshimal 120071 ). 
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5. A lattice of quantized vortices in an atomic BEC 



We now address the issue of a rapidly rotating BEC where many yortices have 
been nucleated and ar ranged into a regular triangular lattice (lAbo-Shaeer 
200ll . ICoddington!l2004l ) . We first present a discussion of the equilibrium prop- 
erties of a rapidly rotating condensate, and then present the basic theoretical 
background for its description. The equilibrium properties of vor tex lattices in 
a trapped BEC have been extensively studied by the JILA group (jSchweikhard 
2004al . ICoddingtonl 120041 ) . We next discuss the collective dynamics of an as- 
sembly of vortices in a trapped BEC. Finally, we discuss an unconventional 
vortex phase which occurs in the presence of an externally applied potential 
created by laser beams. 



5.1. Equilibrium properties 

For very large Q, the rotation of the superffuid mimics a rigid body rotation 
with V x v s = 2fl by forming a vortex lattice. Using the fact that the vorticity 
is given by the form V x v s = k5^- 2 ' (r±)z, we find that the average vorticity 
per unit area is given by V x v s = Kn v z, where n v is the number of vortices per 
unit area. Hence, the density of th e vortices is related to the rotation frequency 
Q as n v = 2Q/k (iFeynmanl Il955l ) . This relation can be used to estimate the 
maximum possible number of vortices in a given area as a function of Q. As 
shown below, the properties of a vortex lattice can be characterized by the 
nearest-neighbor lattice spacing ~ b = [h/mVL] 1 ^ 2 , defined by the area per 
vortex n^ 1 = irb 2 , and by the radius of each vortex core ~ £. 
Note that the GP energy functional of Eq. (11) in a rotating frame can be 
rewritten as 



E' 



dr 



2m 



m 



-iV - — O xr|$ 

n 



cff 



^l 2 + -l^l 4 



(14) 



where V e R = m{u\ — Q 2 )r 2 /2 + muj 2 z 2 /2 is the effective trapping potential 
combined with the centrifugal potential; the rotation effectively softens the 
radial potential and vanishes at Q = u±. Because the first term in Eq. (14) 
reads ft 2 (V|\l/|) 2 /2m-|-m(v s — f2 x r) 2 |\I/| 2 /2, it can be neglected in the TF limit 
and the rigid-body rotation limit v s = f2 x r. Then, the TF radius is given by 
R±(Q) = R±/[l — (f2/co>_|_) 2 ] 3//1 ° with R± for nonrotating condensate and an 
aspect ratio of A r b = R±(Q)/R Z = A/[l — ( fl/u 1 1 2 ] 1 / 2 . Thus, measur ing A r b 



will give the rotation rate of the condensate (lRamanll200ll . lHalj anll200ll ) . Also, 
in the high rotation limit Q — ► u±, the condensate flattens out and reaches an 
i nteresting quasi-2 D regime; current experiments have reached Q/uj± ~ 0.995 
(ICoddingtonll2004h . 



Abrikosov triangular lattice of quantized vortices have been observed experi- 
mentally, as shown in Fig. 2. Direct imaging allows a detailed investigation of 
the nature of a vortex lattice in an inhomogeneous superfluid. Here, we sum- 
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marize the salient results and theoretical interpretations of these observations. 



5.1.1. Lattice inhomoqeneity 



Experimental observations (lAbo-Shaeerll200ll . lEngelsll2002l ) and numerical sim- 
ulation of the 3D GP equation (lFederll2001bl ) have revealed that for a finite-size 
trapped BEC, the vortex density in a lattice is lower than the rigid-body esti- 
mate and the lattice is remarkably regular. S heehy and Rad zihovsky explained 
these points analytically in the TF limit (ISheehvl l2004al lbh : they derived a 
small, radial-position-dependent, inhomogeneity-induced correction term to 
the vortex density as 



n v (r) 
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e 1 k 



(15) 



This result indicates that the vortex density is always lower than the rigid- 
body estimate of the first term in Eq. (15). Also, the vortex density is higher in 
regions where the condensate density is most uniform, that is, the central part 
of a harmonically trapped gas. However, the position-dependent correction is 
small (n v changes less than a few % over a region in which the atom density 
varies by 35%), which seemingly causes regularity of the lattice. These res ults 
have been confirmed by a detailed experimental study (jCoddingtonl 120041 ) . It 
should be also noted that inhomogeneity in the are a density of vortices can 
also be derived in the limit of the lowest Landau level (jWatanabdl2004l . ICooper 
2004] . lAftaliorJlioOol lBavmll2007h . as explained below. 



5.1.2. Attainment of the lowest Landau level regime 

Note that the first term in Eq. (14) can be identified as the Hamiltonian 
Hl = (—iftV — eA/c) 2 /2m of a charge — e particle moving in the xy plane 
under a magnetic field Bz with a vector potential A = (mc/e)Q x r. If the 
interaction is neglected (g = 0), the eigenvalues of the Hamiltonian of Eq. (14) 
forms Landau levels as € n>m /h = u± + ti(uj± + Q) + m(cj_|_ — Q), where n is 
the Landau level index and m labels the degenerate states within a Landau 
level. The lowest energy states of two adjacent Landau levels are separated 
by h(u± + Q), whereas the distance between two adjacent states in a given 
Landau level is h(u± — Q); when Q = u±, all states in a given Landau level 
are degenerate. Physically, this corresponds to the case where the centrifugal 
force exactly balances the trapping force in the x-y plane, and only the Coriolis 
force remains. The system is then invariant under translation and hence has 
macroscopic degeneracy. This formal analogy has led to the predic tion that 
quantum Ha l l-like properties would emerge in rapidly rotating BECs (IHo 2001 



Vieferdl2000L ICooDejl200lilSinovdl200i lRegnauld l200^ I2004L ICazaMa 
Chanel 120051 . lRezavill2005l . iMorrisI 120061 



2005 



Interaction effects mix different (m, n) states. Because the density n of the 
system drops as Q — > u^, the interaction energy ~ gn can become small com- 
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pared with the Landau level separation Ihio^. In this limit, particles should 
condense into the lowest Landau level (LLL) with n — 0. Then, the system 
enters the "mean field" quantum Hall regime, where the wave function can be 



described by only the LLL orbitals with the form = A Y[Az — Zj)e 



where z = x + iy, Zj is the positions of vortices (zeros), and A is a normaliza- 
tion constant. The minimization of Eq. (14) using the ansatz is a use- 
ful theoretical pr e scription to t a ckle the properties of rapidly rotating BECs 
flWatanabd I2004L ICooperl I2004L lAftalionl liooil . ISoninl l2005bl . lAftalionl 120061 . 
Cozzinill2006h . 

Schwaikhard et al. created rapidly r otating BECs by spinning condensates 
to > 0.99 (jSchweikhardl l2004al ) . When the condensate enters the LLL 

regime, characteristic equilibrium properties appear, as described below. 



(i) Global structure 

In the LLL limit, the radial condensate density profile with uniform vortex dis- 
tribution has been predicted t o change from a parabolic TF profile to a Gaus- 
sian profile as Q is increased (lHoll200ll ). However, no signs of such crossover 
have been found in experiments; even when the dynamics were restricted to 
th e LLL, the density profile remained a parabolic TF profile as Q —>■ u)j_ 



(jSchweikhardl l2004al . ICoddingtonl 120041 ) . This result can be seen qualitatively 
from the energy m i nimization under the LLL limit w ith nonuniform vortex 
density (IWatanabd 12004! . ICooperl 12004! . lAftalionl l2005h . As long as the total 
number of vortices is much larger than unity (N v ^> 1), the energy in the LLL 
regime is given by 



E' = QN+ dr 



(u;_l - Q)—n(r) 



bg 



-n r 



'ho 



(16) 



plus terms involving the trapping potential in the ^-direction. Here, n(r) = 
(I^lllI 2 ) is the coarse-grain averaged density profile in order to smooth the 
rapid variations at the vortex cores. Then, the interaction parameter g is 
renormalized to bg, where b = (|^lll| 4 )/(|^lll| 2 ) 2 is the Abrikosov param- 
eter. The energy (16) is then minimized by the TF profile, n(r) = \p — Q — 

(^-ny/aU/bg. 

Since the energy (16) depends only on the smoothed density, the vortices 
adjust their locations so that the smoothed density becomes an inverted 
parabola. In the LLL regime, the relation between the condensate density n(r) 
and the mean vortex density, n v (r) = Y. j S(t — Tj) is given by 4 _1 V 2 lnn(r) = 
~ a ho ~^ 7Tn v( r ) (IH0II2OOII . lAftalionll2005l ) . If the density profile is Gaussian, the 
vortex density is constant. However, for a TF profile, 



n v (r) 



2Q 



■Km 



\2 ' 



(17) 



This result is similar to Eq. (15) in the low rotation regime, where the coef- 
ficient of the second term is different. Since the second term is smaller than 
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the first, as ~ a1 Q / R± — N~ x , the density of the vortex lattice is basically 
uniform, consistent with the argument in Sec. 5.1.1. Turning the argument 
around, very small distortions of the vortex lattice from perfect triangular can 
result in large changes in the global density distribution such that the TF 
form is energetically favored rather than the Gaussian. 



(ii) Vortex core structure and fractional area 
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Fig. 8. Fraction of the condensate surface area occupied by vortex cores A (see text) 
versus the inverse of the LLL parameter rj~^ L , measured after condensate expansion. 
The data clearly show a saturation of A as Q/uj± — * 1. The dashed line represents 
the prediction for the pre-expansion value in the case of a low rotation rate. The 
dotte d line shows the res ults for a saturated value of A in the LLL limit. (Taken 
from ( Schweikhard 2004al ) . Reprinted with permission from APS.) 



Another interesting characteristic of the LLL is that the vortex core is of the 
same size as the distance b = {h/mVL) 1 / 2 between adjacent vortices. The ra- 
dius of a single vortex core is of order £ = 1/ y&tma, so that a vortex core 
would begin to overlap the next at £ ~ b. This gives an upper-critical rotation 
frequency Q c i ~ 8imah/m ~ 10 3 — 10 5 rad/sec, which is an experimentally 
accessible rate. However, there is no phase transition associated with vor- 
tex cores overlapping in a rotating condensate. Rather, vortex cores begin to 
shrink as the intervortex spacing becomes comparable to the healin g length 
and eventually the core radius sc ales with the intervortex spacing ({Fischer 
20031 . lBavmll2004bl IWatanabell27)06h . 

Figure 8 shows the measured fractional area, defined as A = r 2 /b 2 = n v 7ir 2 , as 
a func t ion of the inverse of the LLL parameter T L ll = n/2Ml (jSchweikhard 



2004al . ICoddingtonl 120041 ) . The linear rise of A at small Q occurs because 
the core size remains constant, while n v increases linear ly with Q. Explicitly 
the core radius was estimated numerically as r v = 1.94c; (jSchweikhardll2004al ). 
and n v = mQ/irh by neglecting the effect of inhomogeneity. These values yield 
A = 1.34Flll) shown by the dashed line in Fig. 8, where n = 0.7n pea k and 
/i = gn pea k were used for the estimation. The flattening of A with increasing 
Q is a consequence of the vortex radius scaling with the intervortex spacing. 
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The upper dotted line shows the upper limit A = 0.225, obtained by using the 
oscillator p-state structure l^coreMI 2 ~ [(r/6) exp(— r 2 /26 2 )] 2 as the profile 
of a vortex core in the LLL limit, where b is regarded as the radius of the 
(cylindrical) Wigner-Seitz cell around a given vortex. The data in Fig. 8 show 
the expected initial linear rise, with the predicted scaling of the core radius 
with intervortex spacing. More detailed theoretical studies which treat the core 
structure expl i citly obtain exce llent agreement with the experimental results 
flCozzinill2006l . IWatanabel[2006h . 



5.2. Collective dynamics of a vortex lattice 



5.2.1. Vortex lattice dynamics coupled with collective modes 
As shown in Sec. 4.2.2, vortex states undergo interesting responses to excita- 
tion of the transverse quadrupole modes with m z = ±2. Similar studies have 
been made for rapidly rotating BECs. In this case, the dispersion relation 
of the quadrupole modes is given by u±2 



Q? ± fi (ICozzinil 120031 ) . 
which has been measured experimentally (lHalj an! 120011 ). When Q — * u±, we 
have uj + 2 — ► 2u± and oj-2 — ► 0, reflecting the tendency of the system to become 
unstable against quadrupole deformation. Excitation of the quadrupole mode 
for Q — > u± induces large de formations of the condensate and nonequilibruim 
dynamics of vortex lattices (jEngela |2002| ) . Interestingly, when the m z = —2 
mode was excited, a vortex lattice was distorted to form a one- dimensional 
set of closely spaced vortices. This observation was explained by the fact that 
vorti ces should follow the stream line of the background quadrupole velocity 
field ( )Cozzinill2003l . lMuelleril2003l ) . In contrast, excitation of an m z = +2 mode 
dissolved the regular lattice, where the vortex lines were randomly arranged 
in the x-y plane but were still strictly parallel along the z-axis. 
As stated in Sec. 5.1, a centrifugal force distorts the cloud into an extremely 
oblate shape, and thus the rotating cloud approaches the quasi-2D regime. 
Excitation of an axial breathing mode (m z = 0) has been use d to confirm 
the 2D signature of a rapidly rotating BEC (jSchweikhardl l2004al ) . For a BEC 
in the axial TF regime, an axia l breathing fr equency wb = V^^z has been 
predicted in the limit fl/u± — > 1 (jCozzinill2003l ). whereas uj-q = 2u z is expected 
for a noninteracting gas, expected for \i < hu z . Schweikhard et al. observed 
a crossover of cub from y/3cu z to 2uj z with increasing Vt (/i ~ 3fko z ). Also, 
excitation of the scissors mode in a condensate with a vortex lattice induces a 
collective til ting mode o f the vortex array (the lowest-energy Kelvin wave of 
the l attice) (jSmithl 12004 ) . referred to as an anomalous scissors mode (j Chevy 
2003h . 



5.2.2. Transverse oscillation of a vortex lattice: Tkachenko mode 
The dynamics of vortex lattices itself raises many interesting problems. It 
should be possible to propagate collective waves in a transverse direction to 
the vortex lattice in the superfluid, called Tkachenko (TK) modes. For an in- 
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compressible superfluid, the dispersion law is given by c^xk(^) — \JhVt/Amk. 
The TK modes of a vortex lattice in a trapped BEC have been analyzed theo- 



2004, 
2003) 



retica llv (iBayml l2003l. iMizushima 



20041 . iBaksmatvlliooi lBavml l2004a1. ICozzini 



Giffordll2004l . ISoninll2005al Jbl) and observed experimentally (jCoddington 



Experimentally, TK modes have been excited by the selective removal of atoms 
at the center of a condensate with a resonant focused laser beam, or by the 
insertion of a red-detuned optical potential at the center to draw atoms into 
the middle of the condensate. The former method has also been used to create 



long-lived vortex aggregates (jEngeld 120031 ) . In the experiment, the TK modes 
were identified by the sinusoidal displacement of the vortex cores with the ori- 
gin at the center of the condensate; see Fig. 9(A). TK modes can be classified 
by the quantum number (n, m), associated with radial and angular nodes, in 
a presumed quasi-2D geometry. 

To explain the observed frequency of the TK mode k>(„ jm ), the effects of com- 
pressibility should be taken into acc ount. Accord ing to the elastohydrody- 
namic approach developed by Baym (iBayml 120031 ) . the TK frequency is de- 
scribed by the compressional modulus C\ and shear modulus Ci of the vortex 
lattice, included in the elastic energy 



E P 



J dr j2C 1 (V-e) 2 + C 2 



O^l] \ ( dex \ d€y 
dy I \dy dx 



where e(r,t) is the continuum displacement field of the vortices from their 
home positions. In the incompressible TF regime, Ci = —C\ = nhQ/8. Then, 
the upper branch of the energy spectrum follows the dispersion law u 2 ^ = 

4Q 2 + c 2 k 2 with sound velocity c = ^Jgn/m, being the standard inertial mode 
of a rotating fluid and having a gap at k = 0. Conversely, the low frequency 
branch corresponds to the TK mode and has 



m c 2 k 4 

4m 4Q 2 + c 2 k 2 ' 



(19) 



For large k, this reproduces the original TK frequency c^tk, while for small k 
it exhibits the quadratic behavior ~ J h/lQmVtck 2 . The transition between 
k 2 and k dependence occurs at k ~ VL/c > Rj 1 . This suggests that the effects 
of compressibility, characterizing the k 2 dependence, play a crucial role in the 
TK mode. Thus, this regime is distinguished from the usual incompressible TF 
regime as the "soft" TF regime. When the finite compressibil i ty is inc l uded, 
the observed values of ct>(i,o) ar e well explained (IBayml 120031 . ICozzinil [2004J, 
Soninll2005al ). First-principle simulations based on the GP formalism also ag ree 
excellently with the experimental data (jMizushimall2004l . lBaksmatyll2004l ). 
In the LLL limit, corrections to the elastic shear modulus Ci of the vortex 
lattice are impo rtant; with i ncreasing Q, its value eventually reaches C. 



LLL 



(81/80ir A )mc 2 n ( IBayml 120031 ). Using this value, Schweikhard et al. compared 
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Fig. 9. (A) TK mode (n,m) = (1,0) at A = 1.5 x 10 5 and Q = 0.989u; ± , fitted 
by sine fits. (B) Comparison of measured TK mode frequency uno) (solid symbols) 
versus the theoretical value ( Bavml 120031 ) . using the vortex lattice shear modulus 
Cj F in the TF limit (circles) and Q^ LL in the LLL regime (stars). Note that both N 
and Tlll decrease as Cl = Q/lo±_ increases. For r LLL ~ 3 (reached at A = 7.8 x 10 5 
and Q ~ 0.978) the d ata cross over from the TF to the LLL prediction. (Taken from 
(|Schweikhardll2004ah . Reprinted with permission from APS.) 



the measured u>(i,o) with the theoretical prediction, finding a crossover of ^(i,o) 
from the TF results to the LLL results, as shown in Fig. 9(B). However, more 
detailed theoretical analysis in the LLL limit revealed that the value of the 
shear modulus i s estimated as 0.1027mnc 2 , which is a fac tor of 10 larger than 



Cozzinill2006l ). This indicates that 



(81/807r 4 )mnc 2 ( jSinoval 12002 . ISoninl 120051 
even though the equilibrium properties in the experiment are consistent with 
the LLL picture, the data of the TK frequency are still far from the LLL 
limit. One possible explanation of this discrepancy is an underestimate of the 
rotation rate from the aspect ratio of the cloud d ue to the defocu s of the 



imaging camera or the breakdown of the TF theory ( iWatanabd 120071 ) . 



5.3. Vortices in an anharmonic potential 



For a rotating condensate with a frequency Q in a harmonic potential {l/2)muJ\r i 
the centrifugal potential cancels the confinement, thus preventing a BEC from 
rotating at Q beyond u±. This restriction can be avoided by introducing an 
additional quartic potential, so that the combined trapping potential in the 
xy plane becomes V ex (r) = (l/2)mcc; 2 L (r 2 + Ar 4 /^), where the dimensionless 
parameter A characterizes the relative strength of the quartic potential. The 
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The vortex phases in an anharmonic trap are quite different from those in a 
harmonic trap, since it is possible to rotate the system arbitrarily fast. The 
predicted phase diagram of the vortex states as a function of interparticle in- 
teraction strength versus rotation rate is shown in Fig. 10. For small Q, the 
equilibrium state is the usual vortex lattice state. As fl increases, the vortices 
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Fig. 10. Phase diagram of the vortex state with an additional quartic potential 
for A = 1/2 as a function of interparticle interaction strength (g = 4-7rn z a) versus 
rotation rate (Q c = Q). The dashed curve denotes the onset of a central density 
hole (VLH) in the uniform v ortex lattice state (VL), obtained by TF analysis as 
Ul = 1 + 2^/A(3v / Ag/2vr) 1 / 3 (rFetterjl2005l b The dashed-dotted curve and the solid 
points (•) joined by solid lines show the phase boundary f2 c between the annular 
condensate with a circular array of vortices (AA) and the gia nt vortex (GV) state, 
which are determined by two different analytical methods (Fu 200 6|). The open 
triangles (A) are the values of J7 C determined by the GP solution ( Fett er 2005) and 
the o pen squares (□) are the results using an improved variational approach (|Kim 
20051 ) . The filled square (■) with err or bars gives the a pproximate bounds on f2 c 
determined numerically for g = 125 (jKasamatsul 2002a). F or a we akly interacting 
l imit g ~ 10, a much richer structure was revealed (Jackson 2004al lbl). (Taken from 
(jF\3 l2006h . Reprinted with permission from APS.) 



begin to merge in the central region and the centrifugal force pushes the parti- 
cles towards the edge of the trap. This results in a new vortex state consisting 
of a uniform lattice (multiple circular arrays of vortices) with a central density 
hole. The central hole becomes larger with increasing Q, and the condensate 
forms an annular structure with a single circular array of vortices. A further 
increase of Q st abilizes a gia nt vortex, where all vortices are concentrated in 
the single hole (jFischerl 120031 ) . 

A combined harmonic-plus-quartic potential was formulated by t he ENS grou p 
by superimposing a blue detuned laser with the Gaussian profile ( ]Bretinll2003l ). 
Since the waist w of the beam propagating along the z-axis is larger than the 
condensate radius, the potential created by the laser Uq exp(— 2r 2 /w 2 ) can 
be written as U(r) ~ Uq(1 — 2r 2 /w 2 + 2r 4: /w A ). The second term leads to a 
reduction of the transverse trapping frequency u± and the third term provides 
the desired quartic confinement, giving uj±/2ti = 65 Hz and A ~ 10~ 3 for 
V ex (r). Figure 11 shows experime ntal images o f the condensate density as the 
rotation frequency £7 is increased ( Bret in 20031 ). For < lu_l, the vortex lattice 
is clearly visible. When Q > u±, however, the vortices becomes gradually 
difficult to observe and the images become less clear for Q = 1.05o;j_ (= 2tt x 68 
Hz), which suggests a transition into a new vortex phase. 
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Fig. 11. (a) Density profiles of a rapidly rotating condensate in a quadratic plus 
quartic potential for various stirring frequencies £l/2ir. For these data lv±/2it = 65 
Hz. (b) Ground state structure obtained b y numerical s imulations with a parameter 
corresponding to the experiment of Ref. ( Bretinl 20031 ). The rotation frequency is 
n/27r = 60, 64, 66, 70.6, 73 Hz (respectively, Sl/u± = 0.92, 0.98, 1.01, 1.08, 1.11) 
from left to right, where Q,/2it = 70.6 Hz corresponds to fi^. The first two rows show 
3D views of the vortex lattice as isosurfaces of low atomic density. In the bottom 
row, the d ensity distribu tion is integrated along the z-axis. (Taken from (IBretin 
2003h and (jDanailalboolih . Reprinted with permission from APS.) 



Despite the visibility of the cores, the angular momentum of the condensate 
monotonically increased, confirmed by the measurement of -Rtf and surface 
wave spectroscopy. Hence, the most plausible explanation for this mysterious 
observation is that the vortex lines are still present, but strongly bent when 
Q > u±. This bending may occur due to the finite temperature effect on 
the fragile vortex lattice at a high rotation rate; numerical simulations of the 
3D GP equation show that, when looking for the ground state of the system 
using imaginary time evolution of the GP equation, much longer imaginary 
times were required to reach a well ordered vortex lattice for Q > u± than for 
Q < uj±. Compared with the numerical results shown in the bottom in Fig. 
11, the condensate should still have an ordered visible lattice even for Q > u±. 
To observe the density hole at the center, it is necessary to rotate at a slightly 
faster rate than the upper frequency used in this experiment. 
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5.4- Vortex pinning in an optical lattice 



Rotating BECs combined with an optical lattice are an interesting system, 
which has two competing length scales, vortex separation and the periodicity 
of the optical lattice. The structure of the vortex lattice is strongly depen- 
dent on the externally applied optical lattice. Various vortex phases appear 
depending on the number of vortices per pinning center, i.e., the filling factor. 
The JILA g roup formulated a rotating optical lattice using a rotating mask 
( ITung! 120061 ) . Such a rotating optical lattice provides a periodic pinning po- 
tential which is static in the corresponding rotating frame. The authors have 
observed a structural crossover from a triangular to a square lattice by in- 
creasing the potential amplitude of the optical lattice, as shown in Fig. 12. 
Thes e observations are consistent with theoretical studies ( lReijndersll2004l . |Pu 
20051 ). The rotating optical lattice provides new phenomena in vortex physics 
for rotatin g bosons; such as t he realization of a driven vortex system in a peri- 
odic array (IKasamatsull2006l ) or strongly correlated phases in rotating bosons 
flBhatll2006h 
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Fig. 12. Images of rotating condensates pinned to a co-rotating optical lattice at 
n = 0.866w_l with pinning strength U pin /fi = (a) 0.049, (b) 0.084, (c) 0.143, showing 
the structural crossover of the vortex lattice, (a)-(c) show absorption images of the 
vortex lattices after expansion, (d)-(f) are the Fourier transforms of the images in 



(a)-(c). ko is taken by convention to be the strongest p eak; k , 



fesq, and fctr2 are at 



\-v \-/ •-« — -- o l- i.rl 1 Qsq; — 

60°, 90°, and 120°, respectively, from ko- (Taken from ( Tunell2006l ). Reprinted with 
permission from APS.) 



6. Other topics and future studies 

In this section, we discuss other intriguing problems associated with quantized 
vortices in atomic BECs. Since many theoretical works have predicted novel 
properties of vortices under a variety of situations, it is impossible to refer to 
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all studies in this section. We thus select some important issues that have been 
highlighted in the ultracold atom community or are likely to be investigated 
in experiments in the near future. 

6.1. A vortex in an attractively interacting BEC 

Although we have dealt primarily with vortices in repulsively interacting 
BECs, vortices in attractively interacting BECs characterized by a negative 
s-wave scattering length a < have also received theoretical interest. A ho- 
mogenous BEC with a < is never stable due to self-focusing collapse, but a 
BEC in a confining potential can remain stable as long a s the number o f con- 
densed atoms N lies below a critical value N c ~ aho/H ( Bradley! 1995 ). This 
is because the self-focusing can be balanced by the kinetic energy (quantum 
pressure), which tends to defocus the wave function. 

Since a central vortex state reduces the peak density, it may help stabilization 
of a trapped condensate with a < in the sense that it can contain a larger 
number of atoms by suppressing self-focusing. However, for Q < u± excitation 
of a vortex in a harmonically trapped BEC with a < is prohibited by the 
center-of-ma ss motion, wh ich is the lowest energy state for a given angular 
momentum (I Wilkin! Il998l ). Some authors have proposed that the use of an 
anharmonic confinement can support a stable vortex phase as well as regimes 
of center-of-mass motion (ILundhl 12004 iKavoulakid 12004 lGhoshll2004bl . ICollirJ 
20051 ) . The vortex state in this s ystem can be regarded as a ring bright soliton 
with nonzero winding number (jCarr!l2006l ). showing interesting splitti ng dy- 
nami cs dominated by the dynamical instability of the quadrupole mode (ISaito 
2002h . 



6.2. Vortices in dipolar condensates 



BECs of chromium atoms have recently been created (jGriesmaierl 120051 ) . ex- 
hibiting a larger magnetic-dipole moment (fid = 6/xb; t^B is the Bohr mag- 
neton) than those of typical alkali atoms (fid — This opens the door 
for studying the effect of anisotropic long-range interactions in BECs. The 
interaction potential between two magnetic dipoles /z^e separated by r is 
given by V^r) = (/i /U^/47r)(l — 3 cos 2 9)/r 3 , where fi is the vacuum mag- 
netic permeability and e • f = cos 9. Such dipole-dipole interactions con- 
tribute to the GP equation as a nonlocal mean-field potential as ihdip/dt = 
[-h 2 V 2 /2m + V cx + g\i)\ 2 + J dr'V d d(r' - r)|^(r')| 2 ]^- Since the scattering 
length can be tuned to zero by a Feshbach resonance technique, we ca n obtain 
novel quantum ferrofluids dominated by the dipole-dipole interaction (ILahaye 
2007h . 

The principal effect of Vdd on the equilibrium properties of a condensate is to 
cause distortion of its aspect ratio so that it is elongated along the direction 
of the dipoles. This feature affects the stability of a vortex in a dipolar BEC; 
the thermodynamic critical rotation frequency Q c decreases for a condensate 
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in a pancake-shaped trap (u± < u z ), while it increases a ciga r-shaped tra p 



(u± < u> z ), compared to that of a conventional BEC (Sec. 3.1) (jO'Delll 120071 ). 
Interestingly, the critical frequency Q c can become larger than the onset of 
the dynamical instability of a rotating condensate (see Sec. 3.2.2). This is 
an intriguing regime wher e a rotating dipolar BEC is dynamically unstable 



but vortices will not enter (jBijnenll20071 ). Numerical simulations show that the 
structure of vortices has a craterlike shape for e || z and has an elliptical shape 
for e J_ z (|Pulll999l ). Rapidly rotating dipolar BECs possess a rich variety of 
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6.3. Melting state of vortex lattices: beyond the LLL regime 



At sufficiently hi gh rotation r ates, a vortex lattice should melt via quan- 
tum fluctuations (jSinoval 120021 ) and the system should then begin to enter 
a strongly-correlated vortex liquid phase. Exact diagonalization studies for 
small number of bosons have revealed that the ground states exhibit strong 
analogy with the physics of electronic fractional qu a ntum Hall states (jViefers 
2000l . ICooperl l200ll . IPareded l200l iRegnaultl 120031 . I2004J . iGhoshl l2004ah : for 
specific filling factors v = N/N v , i.e., the ratio of the total number and the 
vortex number, the ground state possesses incompressibility characterized by 
the energy gap. For example, at angular momentum L z = N(N — 1), where 
N v = 2N, the exact ground state is an iV-particle fully symmetric Laughlin 
wave function adopted for bosons: \P(ri, r 2 , • • • , r N ) 



-z k fe-^ r t/ 2a - 



with Zj = Xj + ijjj. 

The conditions for the formation of these states have been expressed as v < 
(9(1). It can be seen that an unrealistically high rotation is necess ary to satisfy 
the co ndition; observed filling factors are always greater than 100 (jSchweikhard 
2004ah . which are still deeply within the mean-field GP regime. To overcome 
this difficulty, insertion of a ID optical lattice along the z-direct ion has been 
proposed to enha nce the quantum fluctuations of the vortices (IMartikainen 
2003 . Snoek! 2006 ). Then, the optical lattice divides the condensate into pan- 
cake fractions coupled by a tunneling process between near neighbors and N 
in a single pancake is greatly reduced. Achieving this regime experimentally 
remains an important challenge. 



6.4- Spontaneous vortex generation associated with phase transitions 



In 2D systems with continuous symmetry, true long-range order is destroyed by 
thermal fluctuations at any finite temperature. For 2D Bose systems, a quasi- 
condensate can be formed with a correlation decaying algebraically in space, 
where superfluidity is still expected below a certain critical temperature. This 
2D phase transition is closely connected with the emergence of thermally acti- 
vated vortex-antivortex pairs, known as the Ber ezinskii-Kosterlit z -Thouless 
(BKT) phase transition occurring at T = T B kt (IBerezinskiil Il972l . iKosterlitz 
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19731 ) . For T < T B ktj isolated free vortices are absent; vortices always exist 
only in the form of bound pairs, formed by two vortices with opposite circula- 
tions. The contribution of these vortex pairs to the decay of the correlation is 
negligible, and the algebraic decay is dominated by phonons. For T > Tbkt, 
the free vortices form a disordered gas of phase defects and give rise to an 
exponential decay of the correlation. 
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In the ENS experiment, a ID optical lattice was applied to an elongated con- 
densat e, splitting the 3 D condensate into an array of independent quasi-2D 



BECs (jHadzibabidl2006l ) . The interference technique revealed the temperature 



dependence of an expo nent of the first-ord er correlation function of the fiuctu 



ating 2D bosonic field (iPolkovnikovl 120061 ) . A universal jump in the superffuid 



density characteristic of the BKT transition was identified by observing the 
sudden change of the exponent, where the finite size effect causes a finite- 
width crossover rather than a sharp transition. Surprisingly, the microscopic 
origin of this transition, i.e., whether or not it is a BKT type transition, was 
directly clarified from the image of the interference of the two 2D condensates. 
If isolated free vortices are present in either of two condensates, the interfer- 
ence fringes exhibit dislocations. Such a dislocation has been observed in the 
high-T regio n of the cros sover, supported by the theory using classical field 



simulations (jSimulall2006l ). 



In contrast, the JILA group inserted a 2D optical lattice into a con densate 



to create 2D bosonic Josephson junction arrays (jSchweikhardl 120071 ). Each 
condensate was localized at a site j. Each had an individual phase 6j and was 
separated by a potential barrier from the nearest neighbors. This system can 
be mapped to the XY model, H = —JJ2(j,j>) cos(9j — 9f), where J denotes 
the tunneling coup ling and the sum is restricted within nearest neighbors 
(ITrombettonil 120051 ) . The XY model is expected to exhibit a BKT transition 
at Tbkt — J from free-energy considerations. Direct imaging of vortex cores 
and the systematic determination of J revealed evidence for a gradual increase 
in the number of isolated free vortices at T > J, consistent with the BKT 
crossover picture. 

A related work is vortex formation by merging thr ee uncorrelate d BECs that 
are initially separated by a triple well potential (jSchererl 120071 ) . Depending 
on the relative phases between the condensates and merging rate, vortices 
formed stochastically without applying rotation. This situation is useful to 
clarify the mechani sm of sponta n eous yortex generation through the Kibble- 
Zurek mechanism ( Kibble] 1976 , Zurek 1985 ) during rapid phase transition 
flLeggettlll998l . lKasamatsuil2002bh . 



6.5. Skyrmions in multi- component BECs 



Another important issue in vortex physics is to elucidate the vortex phases in 
multicomponent (spinor) BECs. Multicomponent order parameters allow the 
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formation of various unconventional topological defects with complex prop- 
erties that arise from interactions between different order parameter compo- 
nents. Since it is possible to load and cool atoms in more than one hyperfme 
spin state or more than one atomic element in the same trap, multicomponent 
condensates can be realized experimentally. Such systems offer an ideal testing 
ground for the study of unconventional topological defects; similar structures 
appear in other condensed matter systems such as superfluid 3 He and uncon- 
ventional superconductors, and theories in high-energy physics and cosmology. 
A few experimental w orks have investigated the properties o f composite vor- 
tices in spinor BECs ( ILeanhardtl 12003 . ISchweikhardl l2004bj . ISadlerl l2006h . A 
review of this topic is presented in Ref. (IKasamatsul l2005bl ) and references 
therein. 



6.6. Vortices in Fermion condensates 



Quantum degenerate Fermi gases provide a remarkable opportunity to study 
strongly interacting fermions. In contrast to other Fermi systems, such as 
superconductors, neutron stars, or the quark-gluon plasma, these gases have 
low densities and their interactions can be precisely controlled over a wide 
range by using a Feshbach resonance technique. For small and negative values 
of the scattering length a the equation of state approaches the limit of a 
noninteracting Fermion gas, while for small and positive values the system 
behaves as bosons of tightly-bound molecules. Therefore, we can study the 
crossover from a BEC of molecules to a Bardeen-Cooper-Schrieffer (BCS) 
superfluid of loosely-bound Cooper pairs when an external magnetic field is 
varied across a Feshbach resonance. Recent topics in th is rapidly grow ing field 
are presented in the comprehensive review paper Ref. (jGiorginil 120071 ). 
Decisive evidence for fermion superfluidity was obtained from observations 
of long-lived vo rtex lattices in a strongly-interacting rotating Fermion gas 
(jZwierleinl 120051 ) . Rotation was applied to an ultracold gas of 6 Li atoms in 
\F = 1/2, mp = ±1/2), in a similar way as for conventional BECs (see Sec. 
3.2.1), with magnetic fields covering the entire BEC-BCS crossover region. A 
crucial problem was detecting the vortex cores in the BCS limit, because a 
suffic ient density depletion at the vortex core could not be expected (jNygaard 
20031 1 . In the experiment, the visibility of the vortex cores was increased by 
a rapid sweep of the magnetic field from the BCS to the BEC side during 
ballistic expansion of the TOF measurement. These measurements strongly 
support the existence of vortices before the expansion even on the BCS side 
of the resonance. 

The most striking aspect of this experiment is that it opens up the possibil- 
ity of studying vortex physics in a strongly-coupled fermion superfluid in a 
systematically controlled way. In the strong-coupling limit \a\ — > oo at the 
resonance, called a unitarity limit, the Fermi gas exhibits universal behavior. 
Along this line, several microscopic calculations of the vortex structure, based 
on the BdG formalism, have been carried out (jBulgad 120031 . ISensarmal 120061 . 
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Chienl 120061 . iMachidal 120051 . 120061 ). More detailed studies may provide useful 



predictions of the mysterious vortex properties in high-T c superconductors, 
and eventually those in room-temperature superconductors. 



7. Conclusion 

Quantized vortices in atomic Bose-Einstein condensates constitute an active 
research field, which has drawn the continuous attention of researchers in re- 
lated fields such as superconductors, mesoscopic systems, nonlinear optics, 
atomic nuclei, and cosmology, as well as superfluid helium. Quantized vortices 
in rotating condensates have provided conclusive evidence for superfluidity 
because they are a direct consequence of the existence of a macroscopic wave 
function that describes the superfluid. The direct imaging of the vortex cores 
and lines helps us understand the fundamentals of superfluid dynamics. Es- 
pecially, the inhomogeneous effect caused by a confining potential yields new 
features in both a slowly rotating regime and a rapidly rotating one, not found 
in a bulk superfluid system. The observed phenomena are consistent with the 
prediction of the Gross-Pitaevskii equation without fitting parameters. 
Finally, this volume addresses mainly the topics of quantum turbulence. The 
feasibility of generatin g quantum turbu lence in a trapped BEC is discussed 



by one of the authors (IKobayashil 120071 ) . 
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